The Acoustic Detection of Ultra High Energy 

Neutrinos 



Jonathan David Perkin 

Department of Physics and Astronomy 
The University of Sheffield 




Thesis submitted for the Degree of Doctor of Philosophy in the 

University of Sheffield 

• April 2007 • 



Abstract 



Attempts have been made to parameterise the thermoacoustic emission of parti- 
cle cascades induced by EeV neutrinos interacting in the sea. Understanding the 
characteristic radiation from such an event allows us to predict the pressure pulse 
observed by underwater acoustic sensors distributed in kilometre scale arrays. 
We find that detectors encompassing thousands of cubic kilometres are required, 
with a minimum of 100 hydrophones per kilometre cubed, in order to observe the 
flux of neutrinos predicted by the attenuation of ultra high energy cosmic rays on 
cosmic microwave background photons. The pressure threshold of such an array 
must be in the range 5 ^ 10 mPa and the said detector will have to operate for 
five years or more. Additionally a qualitative analysis of the first acoustic data 
recorded by the Rona hydrophone array off the north-west coast of Scotland is 
reported. 
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Chapter 1 
Introduction 



Chapters [T] and |2] of this thesis summarise the history of our understanding of 
neutrinos and the methods by which they can be detected. Chapters |3l |4] and |5] 
discuss the potential performance of hypothetical large scale neutrino detectors. 
Chapter [6] reports an analysis of data from an underwater acoustic sensor array 
situated in the north west of Scotland in the United Kingdom. Finally, in Chapter 
[Z|a discussion summarising this work and making predictions for future work is 
given. 



1.1 A Brief History of the Neutrino 

Neutrinos are the second most abundant known particles in the Universe after 
cosmological photons. They are produced copiously, for example, by stellar bod- 
ies, like the Sun, throughout their considerable lifetimes and during their spec- 
tacular deaths. There are, on average, three hundred million neutrinos per cubic 
metre of space, billions penetrate our bodies unnoticed every second. However, 
neutrinos remain poorly understood. 

It is because neutrino interactions are mediated only by the weak nuclear and 
gravitational forces that they acquire their mystery. However, it is this very same 
reason that makes them attractive to the astronomer. Whereas charged particles 
are deflected by magnetic fields, photons are absorbed by inter-stellar matter and 
softened by radiation fields, and more massive particles fall deeper into gravi- 
tational wells; the neutrino traverses the Cosmos retaining its energy and direc- 
tionality until it eventually undergoes a collision far away from its place of origin. 
Should such a collision occur on Earth it could yield information from parts of the 
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Universe no other form of light or matter can reach. 

Whilst detection of neutrinos has been achieved, there has been no direct mea- 
surement of their rest mass. In fact, evidence of neutrino mass has only recently 
come to light. The development of our understanding of neutrinos, pertaining to 
neutrino astronomy, is illustrated in the following time-line (adapted from HI): 

1920-1927 Charles Drummond Ellis and James Chadwick et al firmly establish that the 
beta decay spectrum is continuous |I3/||3|1/ contrary to the predictions of a 
two body decay: n ^ p + e~. 

1930 Wolfgang Pauli postulates the existence of a third, neutral particle present 
in beta decay to explain the observed energy spectrum n ^ p + e~ + 

^neutral- 

1933 Enrico Fermi incorporates Pauli's neutral particle into beta decay theory 
and bestows it the name "neutrino" {little neutral one) BSj. The continuous 
beta decay energy spectrum is explained. 

1953 Fred Reines and Clyde Cowan detect neutrinos from the Hanford Nuclear 
Reactor using a delayed-coincidence technique on the reaction: v + p ^ 
n + /S"*" 10. Neutrinos were registered through observation of the photons 
emitted simultaneously by capture of the neutron n and annihilation of the 
positron /3+. 

1957 Bruno Pontecorvo makes the first hypothesis of neutrino oscillation, in this 
instance between neutrino and anti-neutrino states IZl- 

1968 Ray Davis Jr measures the Solar Neutrino Flux in the Homestake Mine. He 
observes a deficit in the number of interactions compared to predictions by 
John Bahcall's Standard Solar Model. This deficit becomes known as the 
"Solar Neutrino Problem" H. 

1976 The tau lepton is discovered at SLAC and through analysis of its decay it 
is concluded that the tau is accompanied by its own unique flavour of neu- 
trino, the tau neutrino 10. 

1986-1987 The Kamiokande experiment confirms the solar neutrino deficit. Then, along 
with 1MB, the first measurements of neutrinos from a supernova (SN) are 



made ||10||, ||TT]|. Neutrino Astronomy is born. 
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1989 LEP constrains the number of light neutrino species to three: electron, muon 
and tau EZlI. 



1998 Super-Kamiokande reports on the flavour oscillation of atmospheric neutri- 
nos, direct evidence for three different neutrino mass states |T3|. 

2000 Twenty four years after the tau lepton is discovered, DONUT reports on the 
first direct observation of the tau neutrino [[T4|. 



2002 The Sudbury Neutrino Observatory detects both neutral and charged cur- 
rent interactions from solar neutrinos, evidence is mounting in support of 
oscillations as a solution to the Solar Neutrino Problem IITSl. 



2002 KAMLAND confirms neutrino oscillations consistent with the observed so- 



lar neutrino deficit, this time with reactor neutrinos | [T6| 



2002 AMANDA produces the most detailed neutrino sky-map ever, with no evi- 
dence (yet) of point sources |[T7|. 



2006-7 The cubic kilometre extension of AMANDA, codenamed IceCube is under 
construction at the South pole IfTSl . ANTARES nears completion on the 



Mediterranean sea bed [T9l. Neutrino Astronomy is in full flight. 



1.2 Neutrinos in the Standard Model and Beyond 

Neutrinos belong to a family of particles called "Leptons" (from the Greek word 
leptos meaning small or fine). Leptons exist in three flavours and each forms a 
couplet of one charged and one neutral particle as illustrated in Table [LTI 



e 




T 


q = 


e 








q = 






Table 1.1: The leptons. 

The Standard Model (SM) of particle physics predicts that neutrinos are mass- 
less, spin half, fermions interacting only via the weak nuclear force. This imposes 
the condition that all neutrinos are left-handed (and all anti-neutrinos are right- 
handed), since, in the massless limit it is impossible to boost to a frame of refer- 
ence where a neutrino with left handed helicity (projection of spin in the direction 
of momentum) has a right handed chirality (the sign of the helicity). 
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The SM has been constructed from Quantum Field Theories (QFTs) that ex- 
plain the unification of the electromagnetic and weak nuclear forces into the 
electro-weak force via Quantum Electro Dynamics (QED) and the interpretation 
of strong nuclear interactions via Quantum Chromo Dynamics (QCD). Within 
these QFTs there is no gravitational component, so, although the SM describes to 
great precision a wealth of experimental data, it remains for this reason, amongst 
many others, an incomplete theory. Furthermore, as has already been intimated, 
it does not contain a neutrino mass term. 

Although it was not realised at the time, the Solar Neutrino Problem was ev- 
idence for neutrino flavour oscillations. A free neutrino propagating through 
space exists in a superposition of fundamental mass eigenstates Vi, V2 and V3 the 
mixing of which gives rise to the three neutrino flavours Ve, and v^. The neu- 
trino flavours are in fact manifestations of weak nuclear eigenstates. This can be 
illustrated in the form of a matrix equation: 



(v. \ 



Ve 
V 



\Vr ) 



( 



Ue3 ^ 
Ur3 ) 



V2 
\V3j 



(1.1) 



Uel Ue2 

where U is the 3x3 unitary Maki-Nakagawa-Sakata (MNS) matrix |20||. This is 
analogous to mixing in the quark sector as described by the Cabibbo Kobayashi 
Maskawa (CKM) matrix ||2T|. The unitary matrix U can be decomposed into three 
rotations: 



U 



/l 
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Cl3 
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(1.2) 



sinOjj and Cj-y 



cosOij are employed. This decom- 
position allows us to describe three-flavour neutrino mixing in terms of six pa- 
rameters: three mixing angles 6i2> ^23 and 613; a complex, CP violating phase Sep', 

The 



mi 



mi 



and, two mass-squared differences Am^2 = yn\ — m\ and Ism^^ 
sign of Am23 determines the neutrino mass hierarchy; it is said to be "normal" if 
Am23 > 0, or, "inverted" if Am\.^ < 0. Experiments in neutrino oscillation look 
for an appearance or disappearance in one of the neutrino flavours as a result 
of mixing. ^23 arid Am23 are traditionally known as the "atmospheric" mixing 
parameters describing v^, oscillations; and, 612 and Ami2 are the "solar" 

mixing parameters, concerned with Vg flavour oscillations. In each case 
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a deficit (or disappearance) of the expected flux of the first flavour is seen as a 
result of oscillation into the second kind. 

Neutrino oscillations have implications for the neutrino astronomer. The flux 
of neutrinos from cosmic accelerators is predominantly the result of pion decays 
into muons and neutrinos, essentially producing no Vt, i.e. (pv^ : (pv^^ : (pv^ = 
1:2:0. However, since the length of propagation from source to observer is 
very much greater than the baseline for oscillation, it is expected that by the time 
extra-terrestrial neutrinos arrive at the Earth they will be evenly mixed: (pv^ : ^v,, : 
(Py, = 1:1:1. 

1.3 Sources of Ultra High Energy Neutrinos 

Many theoretical models have been proposed that predict fluxes of Ultra High 
Energy Neutrinos (UHEi/s). A recent, short review is presented in |[22|: for the 
interested reader an exhaustive review is given by 1123 II . Mechanisms in which 
Cosmic Rays (CRs) are accelerated to UHE (E > 1 EeV, where 1 EeV= 10^^ eV) 
are described as "Bottom-Up" scenarios. The subsequent weak decay of such 
CRs through pp collisions and interaction with ambient, or intergalactic, radia- 
tion fields produces an associated flux of UHE neutrinos. So-called "Top-Down" 
neutrino production can occur through various exotic mechanisms such as: the 
weak decay of GUT scal^ particles produced during early cosmological epochs; 
annihilation of superheavy Dark Matter (DM) particles; and effects due to the ex- 
istence of Quantum Gravity. Furthermore they can be a signature of "beyond the 
standard model" processes such as violation of Lorentz Invariance and enhance- 
ments to the neutrino-nucleon cross section. A brief discussion of Bottom-Up and 
Top-Down processes follows. 

1.3.1 Bottom-Up neutrino production 

Bottom-Up neutrino production relies exclusively on the acceleration of charged 
particles of cosmological origin as an injector to the neutrino flux. For instance, 
once a CR proton has been accelerated to UHE, it can disintegrate into pions 
following a collision with another proton; this subsequently produces neutrinos 
when the pions decay. 

^ above 10^^ GeV Grand Unified Theories (GUT) propose that the couplings of the strong, weak 
and electromagnetic forces converge. 



5 



Chapter 1. Introduction 



The cosmic ray spectrum 



The high energy (E > 1 GeV) CR spectrum, as plotted in Figure 11.11 extends a 
further twelve decades of energy throughout which the intensity reduces by 33 
orders of magnitude. Some 94 years since Victor Hess first discovered CRs their 
origin is still widely debated. The CR flux below E > 1 GeV is dominated by the 
Solar Wind: coronal ejecta that are readily captured in the Earth's magnetic field, 
giving rise to the aurora phenomenon. The interstellar low energy CR flux is thus 
very difficult to determine. Despite being of great interest, low energy CRs will 
not be discussed further. 

Traditionally the CR spectrum has been well described in terms of two main 
features, namely the "knee" and the "ankle". Both are due to a deviation from the 
underlying power law distribution, due to a change in spectral index (see caption 
of Figure [LT]) . The causes of each of them remain a controversial topic. The origin 



of the "knee", seen as a bump in the spectrum at a CR energy around 10^"^ eV, 
remains an unsolved problem, to which some of the proposed answers include: a 
change in CR composition, interactions with a Galactic DM halo, collisions with 
massive neutrinos, and different acceleration mechanisms to produce the CRs. 
The change in flux at the ankle has traditionally been interpreted as a transition 
from Galactic to extragalactic CRs. More recently it has been suggested that the 
contribution of extragalactic CRs is of importance at lower energies, not far above 
10^'' eV ||25l|. Furthermore, extra detail in the spectrum may yet be resolved as 
more experimental data is acquired, not least a potential second knee and perhaps 
a toe. 

The concept behind Bottom-Up UHECR and UHEv generation is illustrated in 
Figure II. 2[ The Waxman-Bahcall (WB) ||26ll bound on the neutrino flux assumes 
that all neutrinos result from CRs accelerated at optically thin (i.e transparent to 
CRs) sites. Subsequently these CRs undergo pp and p7 collisions after they leave 
the acceleration site. This bound is however speculative, as it assumes that all 
the energy of the proton is transferred to the resultant pion following a collision; 
in reality the pions, and consequently the neutrinos, will have lower energy. The 
WB bound can be exceeded by non-accelerator neutrino sources, or if they ex- 
ist, accelerator sites that are opaque to CRs such that they do not appear in the 
observed CR flux. 
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Fluxes of Cosmic Rays 
(1 particle per m^— second) 



Knee 

(1 particle per m^-year) 



(Swordy -U.Chicago) 
...I I I I 




Ankle 

(1 particle per km^— year) 



ml I I I I I I I 



Q 10 11 1? 13 14 IS 16 17 1R 19 70 ?1 

10 10 10 10 10 10 10 10 10 10 10 10 10 



Energy (eV) 



Figure 1.1: Above 1 GeV the cosmic ray spectrum (from ISH) is well de- 
scribed by a broken E^"' power law, with 7 = 2.7 and two distinguishing 
features, caused by a change in spectral index 7, at 10^^ eV, labelled 'knee' 
and at 10^^ eV, labelled 'ankle'. 
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Target photon / nucleon field 



Figure 1.2: Protons are accelerated to UHEs at acceleration sites across the 
Universe. Matter and radiation fields provide target material for pp and p7 
collisions, producing an associated flux of UHE neutrinos. 

Shock acceleration mechanisms 

"Fermi" acceleration is one of the driving forces that underpins the generation of 
UHECRs. The process of Fermi acceleration involves a transference of the bulk 
kinetic energy of a plasma to the individual ions it contains. Ions are acceler- 
ated, over long periods of time, by shock fronts - regions of compression in the 
plasma at the interface of two areas at different pressures. Successive head on 
collisions with a shock front tend to increase the kinetic energy of an ion. This 
mechanism is sufficient to produce the observed power law spectrum of CRs. 
However, one must add some non-linearity, to describe the shock reaction to ion 
acceleration, and thus ensure the energy spectrum does not diverge. Acceleration 
occurs in the relativistic outflows that surround a powerful central engine, such 
as a super-massive Black Hole (BH) in an Active Galactic Nucleus (AGN). Local 
inhomogeneities in the turbulent magnetic field structures confine ions initiating 
a random walk. Whilst the ion velocities are changed by magnetic confinement, 
the kinetic energy is not. It is possible for an ion to receive multiple kicks from 
a single shock if its random walk overtakes the compression front many times, 
eventually leading to relativistic ion velocities. 

A second method of acceleration, resulting from bulk magnetised plasma in- 
homogeneities has been motivated by observational evidence. It is known as 
"coUisionless shock acceleration" because ions are accelerated by the electromag- 
netic field of the plasma rather than by particle-particle collisions. The details 
of magnetic fields throughout the Cosmos remain elusive and continue to be ac- 
tively studied; there is mounting evidence that coUisionless shock acceleration 
can describe the observed CR spectrum from multiple sources. ||27|. 
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Acceleration sites throughout the Cosmos 

Traditionally Supernova Remnants (SNR) have been regarded as the primary 
Galactic acceleration sites for high energy CRs below the knee. The observation 
of radio, optical and X-ray emission accompanying the emission of charged par- 
ticles from such regions has long provided evidence for electron acceleration up 
to multi-GeV energies. Under the assumption that protons and other ions will be 
accelerated in places where electron acceleration occurs, the difficult task of de- 
termining the leptonic and hadronic components of such regions is underway. In 
particular. X-ray images from the CHANDRA satellite support CR acceleration 
in the forward blast wave of a SNR, a thin shell in the outer shock ||28ll. This has 
been determined from various details such as the shape and form of the emission, 
and the observed spectra. 

If SNRs are typical accelerators for Galactic CRs then Active Galactic Nuclei 
are the traditional sites of CR acceleration beyond the Milky Way. The standard 
model of an AGN is comprised of a central super-massive BH surrounded by a 
dusty torus of matter which forms an accretion disk in the equatorial plane. Axial 
relativistic jets of charged particles are emitted as accreted matter falls on to the 
central BH. It is in the AGN jets that shock acceleration takes place. In the case 
where the relativistic outflow is directed towards the Earth an AGN is detected 
distinctly as a gamma ray source and is classified as a Blazar. 

An alternative candidate accelerator is the Gamma Ray Burst (GRB). GRBs 
are the most energetic events observed in the Universe. There are two breeds of 
GRB: short, typically lasting less than two seconds; and long, typically lasting 
tens of seconds. The luminosity of a GRB is of the order 10"*^ J and can exceed 
that of the combined luminosity of all the stars contained within its host galaxy. 
The sudden gamma ray emission is now thought to be the result of a rapid mass 
accretion onto a compact body, resulting from events such as core collapse super- 
novae in stars of a few solar masses, or binary mergers §. The result of such an 
event is a relativistic outflow of charged particles that penetrates the surrounding 
medium, sweeping up and accelerating ions, however the mechanism by which 
these relativistic outflows are produced remains a hot topic. 

A novel, cannonball (CB) model, of GRBs has been formulated from the hy- 
pothesis that following a core collapse SN into a neutron star or BH, an accretion 
disk is formed around the compact body |!29i. A CB is emitted when a large chunk 

^any combination of black hole, neutrino star, white dwarf etc. collisions may be viable 
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of the accretion matter falls abruptly onto the compact object. The GRB photons 
are produced as the CB falls through the ambient radiation of the host SN and 
photons are Compton scattered up to GRB frequencies, see Figure 11.31 Again CR 
acceleration occurs in the shock fronts of the relativistic outflow. 




Figure 1.3: Schematic of the GRB cannonball model: A "plasmoid" is seeded 
by an abrupt collapse of matter from the torus (T) onto the compact body 
(C). A cannonball (CB) is ejected along the relativistic outflow (R), and, as 
it passes through the remnants of the parent supernova (SNR) photons (7) 
Compton Scatter up to GRB energies. 



The GZK mechanism 

A golden channel for neutrino production exists through pion photo-production 
from UHECR protons. This takes place via excitation of the A""" resonance result- 
ing from interaction with 2.7 K Cosmic Microwave Background (CMB) photons, 
as formulated in Equation II. 3[ 

V + 7CMB ^ A+ ^ n + 7r+ (1.3) 

where 

n p + e~ + 

+ Ve + V^ 
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This reaction was postulated first by Kenneth Greisen 1130 II and later, indepen- 
dently by Zatsepin and Kuzmin | l3T| , and has been dubbed the GZK mechanism. 
The threshold energy for pion production by CR protons on photons of energy 
/ X 10-4 eV (the mean energy of black-body radiation at 2.7 K) is ~ 5 x 10^^ eV. 
Some pion photo-production occurs below this threshold because of the high- 
frequency tail of the black-body photon spectrum. The consequence of this pro- 
cess is to limit the source distance of UHECR protons to within ~ 50 MPc, (within 
the extent of our local group of galaxies) the typical attenuation length in the CMB 
photon field at this energy. The flux of GZK neutrinos appears in Figure 12.131 It 
is highly desirable that a neutrino telescope be sensitive to the flux of cosmogenic 
neutrinos produced via the GZK mechanism, since this represents what is essen- 
tially a guaranteed signal, and a smoking gun for EeV scale CR proton accelera- 
tion. Predictions of the flux of GZK neutrinos must satisfy constraints on the dif- 
fuse gamma ray background due to the process: A"*" p + n^, where n'^ 77, 
which competes with the process in Equation ll.3l Such limits are presently set by 
the Energetic Gamma Ray Experiment Telescope (EGRET), a satellite borne detec- 
tor sensitive to gamma rays in the energy range 20 MeV to 30 GeV. Furthermore, 
they are dependent on an assumed initial proton flux. An example of different 
GZK neutrino fluxes is presented in Figure 11.41 Also noteworthy is the photonic 
analogue of this mechanism, whereby CR photons are attenuated by the cosmic 
infrared, microwave and radio backgrounds. 



1.3.2 Top-Down neutrino production 

The apparent lack of suitable candidate sites for particle acceleration up to UHE, 
within our neighbourhood, has motivated the investigation of possible non- accelerator 
sources. Here, UHE neutrinos constitute some of the decay products of arbi- 
trary, super-massive, X particles. X particle production can occur through a 
number of mechanisms, including the decay of topological defects (a sponta- 
neous break in symmetry^ during a transition in phase); metastable super heavy 
(i.e. m > 10^^ eV) relic particles left over from the Big Bang; or very massive 
Dark Matter (DM) particles. Furthermore, extensions beyond SM physics, such 
as extra-dimensional regimes, can provide scenarios for the production of UHE 
neutrinos. 



^an intrinsic property that renders an object invariant under certain transformations 
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Figure 1.4: Cosmogenic neutrino fluxes (from Il32ll ) as a function of energy. 
The thick solid line (max) is for an initial proton flux (x 1/E, assuming that 
the EGRET flux is entirely due to pion-photoproduction. The thin solid line 
(0.2) shows the neutrino flux when the associated photons contribute only 
up to 20% of the EGRET flux. The dashed line indicates an initial proton flux 
oc The experimental points represent the UHECR flux measured by 

AGASA (crosses) and HiRes (asterisks) - see Section 11.41 for further discus- 
sion. 



X particles from topological defects 

Topological Defects (TD) in states of matter occur as a result of symmetry break- 
ing phase transitions. Terrestrial examples include vortex lines in superfluid he- 
lium, magnetic flux tubes in type II semiconductors and disinclination lines in 
liquid crystals. One usually envisages the formation of TDs as a result of ther- 
mal phase transitions. For instance, a ferromagnet acquires domain structures 
as it is cooled through its Curie point and a symmetry is spontaneously broken 
(see Figure II. 5|) . An extensive review of TD models, including their generation 
and topologies, is presented in Bhattacharjee and Sigl's 2000 Paper "Origin and 
Propagation of Extremely High Energy Cosmic Rays" 1133 II . After the Big Bang it is 
natural to assume that there may have been some temperature through which the 
Universe cooled and underwent a similar process, acquiring domain structures 
||34l . This method of formation, however, appears to contradict the inflationary 
paradigm, which is to dilute the concentration of unwanted TDs by introduc- 
ing an exponential expansion of the Universe through an inflationary phase. It 
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magnetisation 
► 




T > T,„,, G = S0(3) 



T < T,„,, G = S0(2) 



Figure 1.5: As a ferromagnet is cooled through its Curie temperature it 
forms domain regions of aligned magnetic field separated by Bloch walls 
(dotted line). Spontaneous symmetry breaking occurs as one requires the 
introduction of a new degree of freedom: net magnetisation to describe the 
state of the system. Hence there is a transition from the SO(3) to SO(2) group 
symmetry. 



has, though, been realised that TD formation can take place after inflation by way 
of non-thermal phase transitions. Such models have been created that allow for 
abundances of quasi-stable TDs that can exist at present in the Universe 1351 . In 
summary, even if the Universe undergoes a period of inflation during an early 
cosmological epoch, TD formation can form at a later time, due to non-thermal 
phase transitions, thus remaining in concordance with the inflationary paradigm. 

Topological Defects can be described in terms of their Higgs field, the under- 
lying quantum field predicted by theory, with which ordinary matter interacts 
in order to acquire mass (an analogous process is that by which massive objects 
moving through fluids experience drag). Generally, if t] is the vacuum expecta- 
tion value (VEV) of a Higgs field in a broken symmetry phase, then an associated 
TD has core size ^ At the centre of a core the Higgs field vanishes, the 

topologies of core centres are given in Table ll.2[ 



Core Topology 


Topological Defect 


point 


monopole 


line 


cosmic string 


surface 


domain wall 



Table 1.2: Core topologies for different defects. 



Far outside of the core, symmetry is broken and the Higgs fields are in their 
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proper ground states. Hence the 'defect' is a core region of unbroken symme- 
try ("false vacuum") surrounded by broken symmetry regions ("true vacuum"). 
The energy densities of the gauge and Higgs fields within the defect are higher 
than outside and they are stable due to a 'winding' of the Higgs fields around 
the cores. Energy is therefore trapped inside the cores and it is in this way that 
the TDs acquire mass. In general the mass of the TD, m^o, is proportional to Tc 
where Tc is the critical temperature of the defect forming phase transition. For 
a monopole mjQ ^ Tc, for a cosmic string the mass per unit length mju ^ T^ 
and for a domain wall the mass per unit area niju ^ T^. For generic symmetry 
breaking potentials of the Higgs field, Tc ^ tj II33II . Topological Defects can be 
envisaged as trapped quanta of massive gauge and Higgs fields of the underly- 
ing spontaneously broken gauge theory. Sometimes there are quanta of fermion 
fields trapped in the defect cores because of their coupling to the massive gauge 
and Higgs fields, this combination of fields effectively constitutes a massive ob- 
ject contained within the defect - this is the X particle. 

Metastable superheavy relic particles as X particles 

Metastable Superheavy Relic Particles (MSRPs) are an expected bi-product of in- 
flation. In certain models MSRPs can exist in the present with sufficient abun- 
dance as to act as a source of non-thermal, superheavy dark matter p3|. A well 
cited description of CR production from relics of inflation is given in |[36|. One 
has to overcome the problem of producing a particle that has a lifetime that is 
both finite and long enough for it to survive to the present cosmological epoch. 
This paves the way for the introduction of exotic physics such as quantum grav- 
ity and instantons. Any model put forward must adhere to present day limits 
set by experiment. It suffices to say that until there are higher statistics for the 
most energetic events, there will exist an abundance of theories proposing MSRP 
progenitors to the UHECR and UHE neutrino fluxes. 

1.4 Observations of Trans-GZK CRs and Limits on 
the UHECR Flux 

As stated in Section [1.3.11 the threshold energy for the GZK mechanism is ~ 5 x 
10^^ eV. To date, AGASA, AUGER and HiRes have all reported the existence of 
particles above this threshold. The UHECR flux they have measured is shown 
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in Figure |1.6[ The appearance of CRs above the GZK threshold suggests that the 
origin of such CRs is within our local group of galaxies despite there being no 
known sources, or more specifically acceleration regions contained therein, that 
can produce CRs at such energies. As the number of events recorded at these 
energies slowly increases it should become apparent if they are a component of 
the diffuse CR background or they originate from point sources. 
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Figure 1.6: The flux of UHE Cosmic Rays (from Ii37i ) as measured by 
AGASA, AUGER and HiRes. The highest energy AUGER event (> 140 EeV) 
does not appear in this plot because the core of the shower fell outside the 
surface detector and as such it did not pass all quality cuts for flux calcula- 
tions. 



1.5 Summary 

It has been shown in this chapter how neutrinos come to reach the Earth from var- 
ious cosmological origins, most of which are not fully understood. The discovery 
of neutrino mass and subsequently flavour oscillations means that although one 
expects muon type neutrinos to outnumber electron type neutrinos by two to one, 
with no tau neutrino component; in fact, by the time of their arrival at Earth they 
are evenly mixed into equal fractions. 

The flux of UHE neutrinos in which the Earth bathes is intimately linked to 
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the acceleration and subsequent decay of UHECRs. Wherever proton accelera- 
tion has occurred, one can expect to encounter the production of neutrinos with 
comparable vigour. The electromagnetic and hadronic components of cosmolog- 
ical accelerator outflow are retarded by inter galactic matter and radiation. This 
presents a limit on how far we can observe them through the cosmic molasses. 
Neutrinos however can reach us from the furthest depths of the Universe. 

The UHE neutrino flux may not rely exclusively on particle acceleration. The 
lack of suitable acceleration sites within our vicinity has motivated the idea that 
the highest energy CR events that have been observed resulted from the decay of 
very massive entities known only as "X-particles" 
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Neutrino Detection Methods 



2.1 Introduction 

The detection medium for Ray Davis Jr's pioneering Homestake experiment was 
a large tank of tetrachloroethene ||8|. This chemical, essentially a cleaning fluid, is 
sensitive to neutrino capture through the reaction Vgoiar + izCl^'^ e~ + igAr"^"^. 
Neutrinos are emitted by the Sun through the conversion of hydrogen to helium 
and in the decay of Be'^, B^, N^^ and O^^. Amazingly it was by measuring the 
number of individual Argon atoms present after flushing out the detector that it 
was possible to quantify the neutrino flux of solar neutrinos through the appara- 
tus. 

Today, there are numerous experiments across the world either under con- 
struction or in operation, dedicated to the detection of astrophysical neutrinos. 
The principal method by which they are registered is through Cerenkov emis- 
sion. This can be at optical wavelengths via their muon daughters, or at radio 
frequencies via the Askaryan Effect p8|. Large, natural bodies of transparenlBJ di- 
electric media such as Mediterranean Sea water, Antarctic ice, subterranean salt 
domes and even the lunar regolith can serve as natural neutrino calorimeters. 



2.2 Deep Inelastic Scattering 

Figure IZTI illustrates how high energy neutrinos I, interact with a target nucleon 
N, by deep inelastic scattering (DIS) on its constituent quarks. The basic interac- 

^i.e. transparent to the secondary radiation used to observe the neutrino interaction 
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tion is given by Equation 12. 11 

l + N^l' + X (2.1) 

where N is a nucleon, I' is the outgoing lepton and X is one or more excited 
hadrons. Interactions can be neutral current (NC), as mediated by the weak- 
neutral Z boson, or charged current (CC) as mediated by the weak-charged 
bosons. The contributions of these two components to the total neutrino-nucleon 
cross section are illustrated in Figure IATT] in Appendix I A. 1[ 




Figure 2.1: Deep inelastic scattering. I is the incoming neutrino, /' the outgo- 
ing lepton; V is the vector boson (Z or W^) and N is the target nucleon. X is 
the arbitrary excited hadronic final state, k, k' and p are the 4-momenta of the 
incoming neutrino I, the outgoing lepton /' and the nucleon N respectively. 



The square momentum transfer of the interaction, also known as hardness, is 
defined as = —q^ = — (k — k')^. Two dimensionless scaling variables can 
be used to describe the kinematics of DIS, each of which represents a physical 
characteristic of the event. Firstly the Bjorken-x variable, defined as: 

X = .rj^ (2.2) 

2(p • q) 

which is the fraction of the nucleon momentum carried by the participating (va- 
lence or sea, depending on E^) quark. Second is the Bjorken-y variable, defined 
as: 

y = ^ (2.3) 

p • k 

which is the fraction of the lepton energy transferred to the hadronic system. In 
the case of a neutrino DIS then this is the fraction of the neutrino energy that goes 
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into the hadronic cascade, the behaviour of which will be discussed further in 
Chapter |3l 

Calorimetric registration of a neutrino DIS can result through the develop- 
ment of the hadronic cascade and, following CC interactions, through propa- 
gation of the charged lepton. Electrons induce an electromagnetic shower that 
develops coUinearly to the hadronic shower, whereas muons, having a longer in- 
teraction length travel through the detector for several metres before decaying. 
In the case of CC tau neutrino interactions a characteristic "double bang" signal 
can be observed: first a hadronic cascade is initiated by the DIS which is closely 
followed by a displaced secondary cascade induced by the decay of the tau lep- 
ton. 

Particle cascades tend to excite atoms in a calorimeter via excitation and ion- 
isation, furthermore there is Cerenkov emission which can be observed in the 
visible or radio frequency band. This is explained in detail for each case below. 



Optical neutrino telescopes are optimally sensitive to the weak interaction of v^i 
neutrinos into muons via the process + N ^ }i~ + X. The energy of the neu- 
trino is shared between the resultant muon, fi, and the hadronic cascade, X, with 
the muon taking between a half and three quarters of the total neutrino energy; 
it then continues along a path that is effectively coUinear (< 0.5° separation for 
Ey > 1 TeV) to the bearing of the incident neutrino. The charge of the traversing 
muon causes the surrounding medium to become polarised. Subsequent depo- 
larisation of the medium results in the emission of Cerenkov photons along the 
relativistic charged particle track. Interference between photons can occur if the 
wavefronts overlap, which is only possible if the muon travels faster than light 
does in the detection medium. Hence the net polarisation of the medium is as- 
symetric. Wavelets interfere constructively to produce the Cerenkov wavefront if 
/3(= v/c) > 1/n, where v is the muon speed and n is the refractive index of the 
medium. The range of the muon in water extends from ^ 5 m at 1 GeV to ~ 7 km 
at lOTeVEU. 

The number of Cerenkov photons produced, N, within a wavelength interval 
d\ and unit distance dx is given by ||2T|: 



V 
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Figure 2.2: Schematic representation of Cerenkov light emission. 



The wavefront formed by coherent emission of Cerenkov photons forms a 
cone of radiation with opening angle dependent on the velocity of the emitting 
particle, given by Equation 12. 5[ Thus if one can reconstruct the angular aperture 
of a Cerenkov light cone, it can be used to infer the velocity of the transient parti- 
cle. One must then assume a rest mass in order to resolve the particle momentum 
and hence it's energy. 

cos» = i^ (2.5) 

Detection of those photons emitted is facilitated by photomultiplier tubes (PMTs), 
devices that convert optical quanta into charge quanta (conceptually the inverse 
of a light bulb). PMTs tend to be optimally sensitive at wavelengths around 
450 nm and one has to be mindful to choose a detector medium in which the 
scattering and absorption of such light is kept to a minimum. 

Traditionally, a neutrino "telescope" views the Cosmos using the Earth as a fil- 
ter to those muons produced in CR induced air showers high in the atmosphere. 
Hence a neutrino telescope built in the northern hemisphere, looking downward 
through the centre of the Earth, will see a view of the southern sky, and vice versa. 
Muon neutrino interactions can occur in the transparent medium surrounding 
the instrumentation or in the Earth's crust below, so long as the muon has suf- 
ficient collision length to allow it to propagate through the detector, preserving 
the directionality of its parent neutrino. The performance of an optical neutrino 
telescope is limited at low energies by the short length of the muon tracks and 
at high energies by the opacity of the Earth. Since the neutrino-nucleon cross 
section increases with energy, as one looks for neutrinos with greater energy the 
Earth ceases to act merely as a filter and instead becomes impenetrable. 
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A summary of existing projects and proposed extensions follows, including, 
at the end of this section, the current limits on the neutrino flux as set by optical 
neutrino telescopes. 



2.3.1 AMANDA 

Undoubtedly the most advanced optical neutrino telescope to date is the Antarc- 
tic Muon And Neutrino Detector Array (AMANDA) II39II . The most recent second- 



phase AMANDA-II [|40| is an array of 667 downward pointing PMTs, occupying a 
cylinder of one kilometre height and 200 m diameter, located 1.5 km deep under 
the 3 km thick, southern polar ice cap. Construction of AMANDA-II was com- 
pleted during the austral summer of 1999-2000 following the initial phases of the 
shallow ice array AMANDA-A in 1993/94, AMANDA-B4 (80 PMTs on 4 strings) 
in 1995/96 and AMANDA-BIO (an additional 216 PMTs on 6 strings) in 1997/98 
as illustrated in Figure |2.3[ Deployment of detector strings is made possible by 
drilling the ice with heated jets of 75 °C water. Operating at a speed of 1 cm s~^, 
it takes approximately three and a half days to bore a 50 ^ 60 cm diameter hole 
to a depth of 2 km. After a bore hole is formed, the strings are dropped into place 
and the ice re-freezes over a period of about 35 hours. 

AMANDA has been in operation for over 5 years in its latest, phase-II, config- 
uration collecting between 7 and 10 neutrinos per diem IITSB . A background rate of 
^ 80 Hz is present due to downward going muons from atmospheric CR interac- 
tions and must be rejected through angular cuts and event reconstruction. Point 
source searches have been performed using the first four years of AMANDA neu- 
trino data and limits on the neutrino flux have been estimated using 800 days of 
data p9|. Further analysis is currently underway. The most significant point 
source excess, located in the direction of the Crab Nebula, has a confidence level 
of only 3.4(r IfTSl . A view of the sky as viewed through AMANDA is plotted in 
Figure I2.4[ Clearly there is a case for a kilometre cubed extension to AMANDA 
if decisive identification of point sources is to be made. IceCube [ jlHl will incor- 
porate the existing AMANDA-II instrumentation and extend to a total of 4800 
active PMTs and furthermore include IceTop - an array of surface based scintilla- 
tor detectors for the purpose of calibration, background rejection and CR studies. 
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Figure 2.3: Schematic of AMANDA phases A (1993/94) and B-10 (1997/98) 
including graphic of an AMANDA optical module (OM) containing PMT. 
From 133. 
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Figure 2.4: The sky as viewed through the AMANDA neutrino telescope 
(from nil) with the coordinates of the crab a. = 19. Ih, 5 = 9.2° highlighted 
by the • marker. 



2.3.2 ANTARES 

The ANTARES (Astronomy with a Neutrino Telescope and Abyssal Environmen- 
tal RESearch) neutrino telescope is currently under construction in the Mediter- 
ranean Sea, off the coast of southern France |jl9l . It consists of twelve 480 m long 
strings each bearing 75 ten inch Hammamatsu PMTs arranged three per storey, 
each enclosed in a pressure resistant glass sphere (PMT+ sphere = optical module 
(OM)), with a declination angle of 45°. This declination helps limit the loss of 
transparency due to the build up of sedimentation on the glass housing to less 
than 1% per year. Each storey has a box of controlling electronics, the Local Con- 
trol Module (LCM) along with the 3 OMs. There are 5 storeys per sector with a 
vertical spacing of 14.5 m and 5 sectors per line, with a horizontal line spacing of 
between 60 m and 75 m. The pointing accuracy is ~ 1.5° at = 100 GeV but 
improves to less than 0.5 ° for energies greater than 1 TeV due to the increased 
length of the muon track and the subsequent increase in Cerenkov photons ||4T|. 

Since March 2005 a Mini test Instrumentation Line with Optical Module (MILOM) 
has been recording data such as environmental conditions and background rates 
due to bio-luminescence and the decay of solute ^*^K. The MILOM has detected 
a prevailing current of lOcms"^ in the west-east direction, a sea temperature 
of between 13.2 °C and 13.3 °C and a sound velocity of about 1545 m s~^. In June 
2005 the baseline rates for light background were measured to be between 90 kHz 
and 140 kHz at a trigger threshold of 0.5 photoelectrons [|4T|. The angular perfor- 
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mance of the ANTARES telescope relies heavily on knowledge of the positioning 
of its OMs. This is provided by a sophisticated acoustic positioning system that 
has been tested on board the MILOM and has so far met the required 10 cm spatial 
resolution. PMT timing resolutions of ~ 0.5 ns have been measured by illuminat- 
ing OMs with pulsed light from LED beacons mounted at the base of the MILOM 
[|4T|. The first fully instrumented line started recording downward going, atmo- 
spheric muons as of March 2nd, 2006. By April 2007 seven fully instrumented 
lines were in operation. It is expected that full deployment of the completed 12 
string detector will be complete before the end of 2007. 

2.3.3 Lake Baikal 

Baikal is the oldest study group still active; in the late 1970s they trail blazed 
their way in tandem with the ill fated DUMAND (Deep Underwater Muon And 
Neutrino Detector) project | |42| providing today's experimentalists with a wealth 
of information. Located at a depth of just over a kilometre, in Lake Baikal, Siberia, 
is the NT200 neutrino telescope ||43|. In operation since 1998, NT200 comprises 
192 OMs arranged in pairs on 8 vertical strings attached to a fixed metallic frame. 
Each OM consists of a 37 cm QUASAR PMT encased in a transparent, pressure 
resistant housing. Readout and control electronics sit in a similar housing and are 
located midway between each storey. All OMs face downward with the exception 
of the eleventh and second OM on each string, which are oriented in the upward 
direction. Successive OM pairs are separated by a 6.25 m vertical gap, whilst the 
8 strings sit on a regular heptagon, of side 18.6 m, with one string at the centre. 
Since 1998 no neutrino events have passed the NT200 cut selection criteria | |43| . 
An upgrade to NT200+, which will incorporate three new strings and increase 
sensitivity to higher energy events, is underway. 

2.3.4 NESTOR 

The Neutrino Extended Submarine Telescope with Oceanographic Research (NESTOR) 
Project works at a mean depth of 4 km in the Ionian Sea, off the southwestern tip 
of the Peloponnesus |44| . Briefly, the apparatus is composed of rigid titanium 
girders arranged in a six point star formation, with one upward and one down- 
ward facing OM at each tip. Here, an OM corresponds to one 15 inch diame- 
ter PMT along with on board DC -DC converter, surrounded by a spherical glass 
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housing that is pressure resistant up to 630 atm. The diameter of each hexago- 
nal floor (star) is 32 m. Aim sphere at the centre of each star houses the local 
electronics. A completed NESTOR tower consists of twelve such elements in a 
vertical arrangement with a separation of 30 m. To date one star element has 
been deployed for testing purposes accumulating over two million triggers in 
the longest period of continuous operation. Proof of principle has been achieved: 
calibration via LED flashers; detection of backgrounds from the decay of K^*^ and 
bio-luminescence; and reconstruction of atmospheric muon tracks all support the 
continued development of the full scale NESTOR detector. 

2.3.5 Kilometre-Cubed optical Cerenkov neutrino telescopes 

Even before the first generation of large scale neutrino telescopes are complete 
a second wave of activity emerges, this time on an almost unimaginable scale. 
The experiments described previously have become test benches for the world's 
largest particle detector experiments - cubic kilometre neutrino telescopes: 

IceCube 

IceCube is the natural successor to AMANDA and, as mentioned in Section |2. 3. 1[ 
will incorporate the existing AMANDA-II apparatus as a means to upgrading to 
an enormous 4800 OMs and encompassing a cubic kilometre of Antarctic ice. 

By the end of summer 2006 IceCube had 9 new strings and 16 of the IceTop 
surface Cerenkov detector stations; a further 858 DOMs had been approved for 
the next stages of fitting. Each of the deployed strings, already acquiring data, 
consists of 60 Digital Optical Modules (DOM) on a 2.5 km cable IITSl . 

KM3NeT 

Funding from the European Union is now in place to undertake the KM3NeT de- 
sign study phase. Building on the experiences of ANTARES, NEMO and NESTOR, 
the KM3NeT collaboration seeks to develop a water based partner to IceCube, at a 
site in the Mediterranean, yet to be determined. Such a telescope, in partnership 
with its austral cousin, would complete a full sky coverage, allowing for some 
overlap in the plane of the Galactic centre. The Galactic centre itself is viewable 
exclusively from the northern hemisphere. 
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NEMO 

The NEutrino Mediterranean Observatory (NEMO) project represents an alterna- 
tive ongoing effort to pilot a cubic kilometre optical Cerenkov neutrino telescope 
in the Mediterranean Sea. Whilst full optimisation of detector geometry has yet 
to be performed, a demonstration model incorporating some of the features of a 
kilometre cubed detector was deployed at a site off the coast of Sicily in December 
2006 1151 



2,3.6 Limits on the neutrino flux from optical Cerenkov neutrino 
telescopes 

Current limits on the neutrino flux from optical Cerenkov based detectors appear 
in Figure 12.51 Some predicted limits are also included based on proposed exten- 
sions to existing experiments and from new projects. Plotted on the same figure 
are some model neutrino fluxes predicted by theory. 



2.4 Radio Cerenkov Detectors 

Radio Cerenkov detectors are sensitive, not to the passage of muons, but to the 
electromagnetic component of neutrino induced particle cascades. The optimum 
process for detection is a charged current Vg interaction. This simultaneously 
initiates an electromagnetic and a hadronic particle shower from the emergent 
electron and hadron respectively. The two showers are in principle superim- 
posed. In 1962 G. A. Askaryan realised that the net charge of the particle cas- 
cade, in a dense medium, will not be neutral - since positrons emitted by pair- 
production will annihilate in flight, leaving a negatively charged excess. Addi- 
tional electro-negativity is provided by the accumulation of delta rays and Comp- 
ton scattered electrons. In his paper [|38|1 Askaryan predicts that at the maximum 
of an E = 10^^ eV shower there is a 10% electron charge excess; one can expect to 
observe a net 20% — 30% negative charge excess. 

Almost forty years after this effect was postulated came experimental con- 
firmation. Formerly, radio emission from cosmic-ray induced particle showers 
was shown to be measured as a result of the geomagnetic separation of electric 
charges. Through the detection of photon induced cascades in 3200 kg of silica 
sand, evidence for the observation of the Askaryan effect was given by corre- 
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Figure 2.5: Constraints on the neutrino flux set by optical neutrino tele- 
scopes. The • markers represent the AMANDAII atmospheric neutrino data 
set, the thick solid line is the 800 days AMANDAII sensitivity and from the 
same paper is an estimate of the Waxman-Bahcall neutrino flux predicted 
from acceleration of the highest energy cosmic rays lfT8| . The latest limits 
from the Lake Baikal experiment are plotted as the thin solid line ||43| . The 
dashed, dotted and dash-dotted curves represent the projected sensitivities 
from one year of ANTARES IMl, KM3NeT [47], and IceCube IH respec- 
tively. 



lation of the radio pulse size with the energy of a shower. Furthermore it was 
shown that the experimental observations were inconsistent with geomagnetic 
charge separation and consistent with Cerenkov emission because of the align- 
ment of the polarisation along the shower axis and not the local geomagnetic dip 

m. 

A brief report of existing radio Cerenkov neutrino telescopes follows, ended 
by a summary of experimental limits they impose on the neutrino flux. 



2.4.1 ANITA 

The balloon borne ANtarctic Impulsive Transient Antenna (ANITA) experiment 
is a recent and innovative apparatus designed for the radio detection of CR neu- 
trinos. An 18.4 day flight has been achieved by a prototype, named "ANITA-lite" 
early in 2004. 1491 The principal ANITA design features a 2n array of anten- 
nas sensitive over the frequency range 230 — 1200 MHz, observing ^ 2Mkm'' 
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of Antarctic ice, where the depth is equivalent to the radio attenuation length, 
from an altitude of 37 km. The ANITA-lite prototype, whilst only consisting of 
two ANITA antennas with 12% of the viewing area of the full scale detector, still 
observed a considerable 0.17 Mkm"^ of ice, producing competitive constraints on 
the UHE neutrino flux that ruled out the remaining "Z-Burst" {vv Z^) flux 
that was not already constrained by previous experiments, as a mechanism for 
the production of the highest energy CRs. Simulations of the full scale detector 
project a sensitivity down to sub-GZK fluxes, thus opening a window to cosmo- 
genic neutrinos and possibly never before seen, extra galactic cosmic accelerators. 
A two month mission successfully launched in December 2006. 



2.4.2 FORTE 



Developed by the Los Alamos and Sandia National Laboratories the Fast On- 
orbit Recording of Transient Events (FORTE) satellite sits in a circular orbit at 
an altitude of 800 km and inclined 70° from the Earth's equator. Launched in 
1997, its on board optical and radio frequency receivers can be exploited for the 
study of hadronic particle cascades in the Greenland ice sheets as well as light- 
ning storms and for its intended use as a nuclear detonation detector. The most 
recent analysis of the FORTE data set contains over 4 million events, recorded 
from September 1997 to December 1999L50J. Only a single event passed through 
each level of cuts, so the calculated flux limits for the experiment assume this to 
be a representation of the background level. Backgrounds are present in the form 
of radio emissions from lightning strikes. Certain artifacts of a lightning event 
can be used to reject backgrounds of this nature such as signal grouping, random 
polarisations and pulse pairs resulting from reflections off the ground. Lightning 
activity over the Greenland is however extremely rare and thus good knowledge 
of the geographic location of a signal can further reduce the background activ- 
ity. Tentative limits have been placed on the flux of astronomical neutrinos in the 
energy range 10^"^ 10^^ GeV. One treats these with caution noting that they ex- 
tend beyond the scale for GUT particle masses and are calculated in the absence 
of neutrino-nucleon cross section data above 30 TeV. 
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2.4.3 GLUE 

As stated at the beginning of this chapter, it is possible to observe neutrino in- 
duced cascades that occur not in sea, ice or salt on Earth but in the lunar regolith; 
this is illustrated schematically in Figure |Z6l The Goldstone Lunar Ultra-high en- 
ergy neutrino Experiment (GLUE) looks for radio signals from this source. This 
technique is facilitated by two large antennas, one 70 m and one 34 m in size, sep- 
arated by a baseline of 22 km. By utilising the two receivers in coincidence, with 
a large separation, one can in principle reduce the anthropogenic radio back- 
ground. On average the background due to random coincidences and thermal 
noise is reported to be 3 x 10""^ Hz, resulting in a > 95% livetime during opera- 
tion El. 



front 
side 



cascade 





neutrino 
enters 
moon 



partial Cherenl<ov cone 



Figure 2.6: UHE cascades in the lunar regolith (from [51]). 



The reported analysis incorporates a livetime of approximately 120 hours which, 
although quite small, employs up to 0.1 Mkm'^ of detector material, which is 
enough to constrain some of the aforementioned Z-Burst and TD models. No 
neutrino candidate events are present in the 120 hours data set IISTl . 



2.4.4 LOFAR/LOPES 

LOFAR, the LOw Frequency ARray, including the LOFAR PrototypE Station (LOPES), 
is a large scale radio telescope, being built for a broad range of astrophysical 
studies. One goal is the detection of CRs and neutrinos through geo-synchrotron 
emission and radio fluorescence in the Earth's atmosphere and the Askaryan ef- 
fect in the terrestrial and lunar regoliths II52II . The operable range of frequencies 
is from 15 240 MHz. The cleanest channel for direct neutrino detection is 
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through Cerenkov emission in the lunar regolith |[53| in the same way as illus- 
trated by Figure 12.61 It is predicted that the efficiency of the technique should 
allow for detection of neutrinos an order of magnitude below predicted GZK 
fluxes, after only 30 days observation at a signal detection threshold of 500 x 
10-26 Wm-^Hz-i ||5l. 



2.4.5 RICE 

RICE is the Radio Ice Cherenkov Experiment, at the South Pole. It was deployed 
in the Antarctic ice in tandem with AMANDA (section 12.3. 1|) and consists of an 



18 channel array of radio receivers distributed about a cube of length 200 m. The 
depth of each receiver varies between 100 — 300 m. Each of the dipole anten- 
nas has a dedicated, in-ice pre-amplifier before the signal transmission along a 
cable to the surface where background noise from local electronics and more 
importantly the AMANDA PMTs is filtered away. Event triggering initially re- 
quires that at least four channels register a peak above a threshold set beyond 
the level of thermal and background fluctuations, or that any channel registers 
a peak above threshold in coincidence with a 30— fold AMANDA-B trigger. A 
horn antenna at the surface provides an active veto against surface generated 
background transient events. Raw trigger rates of 10 Hz are reported 1155 II after 
veto. There is a reported livetime of ^ 90% after discrimination of surface tran- 
sients of anthropogenic origin. Event vertices are reconstructed via the difference 
in arrival times of signals at successive pairs as well as a lattice interpolation algo- 
rithm; a Cerenkov cone of opening angle 57 ° is then superimposed on the vertex. 
A location accuracy of ^ 5 m is reported. Data analysis spanning records from 
1999-2005 has been performed | [56ll yielding no events consistent with a neutrino 
induced cascade. 



2.4.6 SALSA 

The SALt dome Shower Array (SALSA) collaboration represent the greatest in- 
terest in the instrumentation of subterranean salt domes with radio antennae 
for the detection of ultra-high energy, neutrino induced particle cascades. To 
date preliminary investigations of United States sites have been undertaken as 
a means to surveying the suitability of naturally occurring halite deposits for 
use as calorimeters. Indications are that in terms of enclosed mass within one 
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attenuation length radius, salt domes provide a competitive material to ice. Fur- 
thermore, the simple cubic-lattice structure of halite is potentially less causative 
of distortion through depolarisation and rotation of the signal about the plane of 
polarisation. Additionally halite is less birefringent than ice too. Results from a 
Monte Carlo (MC) simulation of an instrumented salt dome have been reported 
II57II . A sensitivity of the order 10 events per year from a minimal GZK flux is 
expected from a 10 x 10 x 10 antenna array on a 200 m grid spacing. Perhaps 
more important is the recent experimental work that has confirmed the genera- 
tion of coherent Askaryan radio pulses, in salt, as a result of a net charge excess, 
synthesised during the development of particle cascades II48B . As discussed at 
the beginning of this section, this work clearly upholds the production of radio 
pulses via the Askaryan effect and refutes the hypothesis of geometric charge 
separation as their source. 

V 

2.4,7 Limits on the neutrino flux from radio Cerenkov 
neutrino telescopes 

Current limits on the neutrino flux from radio Cerenkov based detectors appear 
in Figure 12.71 Some predicted limits are also included based on proposed exten- 
sions to existing experiments and from new projects. Plotted on the same figure 
are some model neutrino fluxes predicted by theory. 

2.5 Acoustic Detection of UHE Neutrinos 

Acoustic detection of high energy particles was first postulated in 1957 1158 II and 
subsequently full theoretical analyses of the signal production mechanism were 
performed (e.g. ||59l). The hadronic particle cascade induced at the interaction 
vertex by an UHE neutrino has enough thermal energy to locally heat the sur- 
rounding medium, causing it to rapidly expand. This produces a pressure pulse, 
measurable on a suitable acoustic receiver. The minimum neutrino energy re- 
quired for this technique lies around 10 EeV (10^^ eV). Whilst the structure of the 
hadronic showers is the same for all three flavours of neutrino, the behaviour of 
the lepton differs. For Vg type neutrinos the lepton energy is deposited as an elec- 
tromagnetic cascade that is effectively detected along with the hadron shower. In 
the case of v^s the resulting muon is virtually undetectable acoustically because 
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Figure 2.7: Constraints on the neutrino flux limits from radio Cerenkov de- 
tectors. The sensitivity of ANITA Lite is shown by the o markers and the 
projected sensitivity of the full ANITA detector by the • markers Il49ll . The 
FORTE sensitivity is shown by the ■ markers EHl. GLUE EH and RICE (SSl 
are indicated by the ▲ and x markers respectively. Again the WB flux of 
neutrinos from IfTSlI is shown. 



the mean free path for catastrophicj^Bremsstrahlung and pair-production is of the 
order of several km. interactions produce a tau lepton which may or may not 
decay within the fiducial volume of a detector, initiating a second cascade. 



2.5.1 The LPM effect 

Before discussing the acoustic detection of neutrinos it is necessary to consider 
the dynamics of such energetic events. At EeV energies quantum interference 
affects the passage of leptons through a medium. The Landau Pomeranchuck 
Migdal (LPM) effect ||60l l ,|6T | describes a suppression of the cross sections for 
pair production and Bremsstrahlung. This effect dominates the development of 
the leptonic component of a neutrino DIS event as its interaction lengths become 
comparable to the interatomic distances of the medium through which it propa- 
gates. An EeV electromagnetic cascade that is typically a few metres in length at 
TeV energies can extend to hundreds of metres because of the LPM effect. The dis- 
cussion of thermoacoustic emission that follows is limited to the hadronic compo- 

^meaning nearly all the particle energy is lost in one collision 
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nent of the neutrino interaction. Any acoustic emission from the electromagnetic 
cascade is assumed to be inaudible due to the extended nature of the source and 
is thus neglected. 

2.5.2 Formation of the acoustic signal 

The thermal energy deposition along the cascade axis occurs via ionisation and 
excitation of the surrounding medium and is instantaneous on the acoustic and 
thermal diffusion timescale. The cascade can be thought of as a series of dis- 
crete regions in which a Gaussian heat deposition occurs, resulting in an instan- 
taneous step in temperature. The corresponding pressure wave is simply the 
second derivative with respect to time of the change in temperature, yielding a 
distinctive bipolar acoustic pulse. Because of the instantaneous development of 
the cascade, the acoustic radiation is emitted in phase such that individual pulses 
interfere. This is analogous to a line of sources emitting Huy gens' wavelets. In the 
far field these sources have undergone Fraunhofer diffraction such that the total 
acoustic radiation is confined to a narrow, so-called "pancake" with an opening 
angle 6 of about 5°. The geometry of the cascade and the acoustic emission is 
illustrated in Figure 12. 8[ 



Figure 2.8: Schematic of the energy deposition. At Ey = 10^^ eV, in water, 
a hadronic shower deposits 99% of its energy in a cylinder of length L = 
20 m and radius R = 20 cm. The energy deposition of the shower can be 
considered as a series of coherent sources lying along the shower axis. 

Following the workings of | |59ll and ||62| the acoustic signal production is anal- 




X 
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ysed. If the energy per unit volume per unit time is given as a function E(r, t), 
then the total neutrino energy is Eq = JyE{r)d^r. The wave equation for the 
pressure pulse produced p is: 



Where for seawater the parameters are: 

c, the speed of sound in water ^ 1500 ms~^, 

jS, the bulk coefficient of thermal expansion ^ 2.0 x lO"'* K~^, 

Cp, the specific heat capacity at constant pressure ^ 3.8 x 10^ Jkg^^K"^, 

coq, the characteristic attenuation frequency ~ 25 GHz. 

It should be noted that coq is actually a function of the frequency of the radia- 
tion (it is, in fact, dependent on the sound attenuation coefficient, which is itself 
frequency dependent) but to keep these calculations simple it is taken as a con- 
stant for the frequency range of interest / = 10 — 100 kHz. One can express the 
instantaneous nature of the heat deposition by E{r',t) = E{r')S{t), where r' is 
the location where the instantaneous heat deposition takes place. The resultant 
pressure wave some other location r as a function of time t is therefore: 

p(r,i) = / E{t')G{T-r',t)d\' (2.7) 

Jv 

G{r,t) is the pressure pulse generated by a point source E = S{r)S{t), taking 
attenuation into account this is: 

G(M) = -j4.^^^.-"-"'>^'P-' (2.8) 



where t = ^/r/ (coqc). 



2.5.3 Laboratory based measurements of thermoacoustic 
emission 

Since the thermoacoustic mechanism for detection of particle cascades was pro- 
posed 1159 II a number of experiments have attempted to measure the effect in the 
laboratory, either through the use of particle accelerator beams or high energy 
light sources. A discussion of two such experiments follows. 



34 



2.5. Acoustic Detection of UHE Neutrinos 



The 1979 Brookhaven-Harvard experiments 

Some twenty-two years after Askaryan's proposal of the thermoacoustic mech- 
anism for high energy particle detection, the first experiments were undertaken 
at the Brookhaven National Laboratory (BNL) and Harvard University cyclotron, 
to record the acoustic signal emitted due to proton beams traversing fluid media. 
The following summary is based on the material reported in 1631. Three experi- 
mental apparatus were used: 

1. The BNL 200 MeV proton linear accelerator (LINAC) produced a beam of 
protons with total bunch energies between 10^^ 10^^ eV. The range of the 
beam in water was ~ 30 cm, with a fixed 4.5 cm diameter. The spill time 
was variable between 3 200 \is. The energy of the beam was not tunable. 

2. The Harvard University cyclotron accelerated protons up to 158 MeV with 
a beam energy that could be tuned down to 10^^ eV (close to threshold for 
the acoustic signal). The beam had a range of 16 cm in water with a min- 
imum 50 las deposition time, which was long in comparison to the sound 
transit time across the diameter of the beam and thus dominated the tem- 
poral structure of the emitted acoustic signal. The beam diameter could be 
varied between 1 ^ 10 cm. 

3. A second fast extracted beam (FEB) from a 28 GeV proton accelerator was 
utilised at BNL. This had a fixed beam energy of lO^*^ eV, with a range of 
20 cm. The diameter of the beam was variable between 5 ^ 20 mm and the 
deposition time was short at 2 ]is. In contrast to the cyclotron experiment, 
where the beam spill time was the dominant effect, the sound transit time 
across the beam diameter dominated the temporal structure of the signal. 

Each of the above apparatus had sufficient dimension to allow for resolution of 
the initial signal and subsequent reflections (much greater than the beam diam- 
eter or cascade length). The primary source of error was due to the uncertainty 
in the measurement of signal amplitudes and beam intensities. It is reported that 
audible clicks could be heard by the unaided ear, but this in itself did not prove 
the thermoacoustic mechanism. It still remained to disprove acoustic emission 
due to molecular dissociation and microbubble formation. 

In the case of a short beam spill time (10 ]xs), where the signal transit time was 
d/c = 4.5 cm/1.5 mm ias~^ = 30ias in water at 20 °C, a clear bipolar signal was 
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registered, this is plotted in Figure 12.91 The separation between this and the sub- 
sequent first reflection displayed the appropriate delay for the given geometry. 
The half width of the pulse was directly proportional to the transit time across 
the diameter of the beam and the time delay between hydrophones at two differ- 
ent locations was consistent with their respective path lengths. A leading com- 
pression was confirmed by placing two hydrophones 180 ° out of phase. When 
the beam spill time was significantly longer (100 las) than the signal transit time 
(30 ]xs) then there was a clear separation between the leading compression and 
the subsequent, final rarefaction. The observed signal (Figure |2.10|) is calculated 
as the convolution of a quasi-instantaneous pulse (Figure |Z9l) , and the beam spill 
intensity (Figure |2.11|) as a function of time. 

The signal period as a function of spill time served to act as a proof of principle 
for the thermoacoustic mechanism and its underlying linear theory. The level 
of agreement was limited to 2a (85%) due to the large errors imposed by the 
uncertainties in the measurements of the beam intensities and pulse heights. 

Repeated tests were performed with different media such as olive oil and CCI4 
to look for deviations from the linear models as a result of excitation of ions 
and /or molecular dissociation. Microbubble formation was also examined by 
using the different media under varying combinations of temperature and pres- 
sure (so as to alter the microbubble diameters). No deviation was observed as a 
result of either ion excitation, molecular dissociation or microbubble formation 
within the 25% experimental uncertainty. Furthermore the spatial dependence of 
signal amplitudes was shown to vary as the reciprocal of the distance away from 
the cascade as anticipated by the model of linear superposition in the far-field 
condition. 

The aforementioned experiment serves as a proof of principle for the thermoa- 
coustic model. There were large inherent uncertainties to the technique used for 
measuring the thermoacoustic emission from each of the accelerator beams used, 
constraining the agreement between the observed data and the predictions of the 
underlying theory only to within 2a. 

Uppsala proton beam and Erlangen laser beam studies 

The most recent attempt to test the predictions of the thermoacoustic mechanism 
||64| have utilised both a high energy (10 400 PeV total bunch energy) proton 
beam and an EeV pulsed infra-red (1064 nm) Nd:YAG laser. The proton beam. 
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Figure 2.9: The pressure pulse observed from a 10 las proton beam spill (from 
\63\ ). The hydrophone producing the lower trace is 180 ° out of phase with 
the upper one. 




Figure 2.10: The pressure pulse observed from a 100 las proton beam spill 
(from [|63|). The hydrophone producing the lower trace is 180 ° out of phase 









f^^^. . - - ^^^m.'^' -m 



Figure 2.11: Time profile of 100 vis beam intensity (from 116310 . 
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being charged, both ionises and excites the target medium but the laser beam, 
being neutral, will only excite target atoms. The density of pure water is maxi- 
mal at a temperature of 4°C; below this the coefficient of thermal expansion (/3 
in Equation I2.6|) is negative. Hence an increase in temperature in pure water be- 
low 4 °C produces a compression and not an expansion. The coefficient for sea 
water is greater than that for pure water and increases with increasing pressure, 
temperature and salinity ||65|. In this experiment the temperature of the target 
water tank was varied between 1 °C and 20.0 °C ±0.1 °C. The proton beam, deliv- 
ered from the 177 MeV cyclotron at the Theodor Svedberg Laboratory in Uppsala, 
Sweden, had a diameter of 1 cm and a spill time of 30 \is. It produced a shower 
with approximately uniform energy deposition in water for 20 cm, ending in a 
Bragg Peak at approximately 22 cm. The Erlangen Physics Institute laser energy 
was adjusted between 0.1 lOEeV with a fixed pulse length of 9 ns and beam 
diameter of a few mm. The energy deposition of the beam, along the beam axis, 
decayed exponentially with an absorption length of 6.0 ±0.1 cm. 

Acoustic pulse data from the laser produced the expected inversion, predicted 
by the thermoacoustic model, at 4.0 °C to within 0.1 °C. The same measurement 
using the proton beam however revealed an apparently temperature indepen- 
dent contribution to the overall amplitude, present at a level of 5% of the signal 
at 15 °C. The results predicted by the thermoacoustic mechanism were obtained 
following subtraction of this artifact from the proton beam data. The source of 
this contribution to the acoustic signal was not verified and was identified as a 
topic for further investigation. A hypothesis was proposed that the extra signal 
was a result of the charge of the proton beam, hence it was not apparent in the 
laser induced pressure waves. 

Whilst the results of the Erlangen-Uppsala tests still leave some unanswered 
questions, they appear to strengthen the support of the thermoacoustic model. It 
still however remains to identify the extra source of sound apparently resulting 
from the charged nature of the proton beam. 

2.5.4 Experimental results from SAUND 

The SAUND collaboration | [66| have published |[67| the first diffuse neutrino flux 
limit based on their work on a United States Navy hydrophone array located in 
the Tongue of the Ocean, a deep basin situated near the Bahamas. They use a 
seven hydrophone subset of the total 52 hydrophone AUTEC array, arranged in 
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a hexagonal pattern at depths between 1570 m and 1600 m. The spacing between 
hydrophones is between 1.50 km and 1.57 km and they are mounted upon verti- 
cally standing booms that are anchored to the seabed. 



1570^ 




Figure 2.12: Schematic view of the hydrophone array used by the SAUND 
collaboration (from [67| ). 

A novel technique was employed by which calibration of the hydrophones 
was achieved by dropping weighted light bulbs into the sea from a stationary 
boat at the surface. Measuring the signal as the bulb imploded under pressure 
gave a rough estimation of the energy sensitivity of the array to a known source 
and allowed for timing calibration. 

In addition to the AUTEC array performance evaluation a computer simula- 
tion was developed to test the sensitivity of two hypothetical arrays. Each was 
comprised of hexagonal lattices of 1.5km long strings, the strings being modelled 
to have continuous pressure sensitivity along their entire length. Array "A" is 
bounded by circle of radius 5km with 500m nearest neighbour spacing. Array 
"B" is bounded by a circle of radius 50km, with 5km spacing between strings. 
The results of this work are illustrated in Fig 12.131 the next phase of this experi- 
ment "SAUND 11" is currently underway. 

2.6 Practical Motivation for the Acoustic Technique 

The introduction hints why an astronomer would look for neutrino light: neutri- 
nos interact only via the weak nuclear and gravitational forces. This means they 
can reach us from the furthest depths of the Cosmos without being deflected, illu- 
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Figure 2.13: Flux limit curves for the AUTEC array (solid line labelled 
SAUND) and the hypothetical arrays (A and B) including the predicted flux 
of several theoretical models: "Z" for Z-Bursts, "TD"for Topological Defects, 
"WB" for the Waxman and Bahcall flux and "GZK" for the GZK theoretical 
flux. Some experimental limits (FORTE, GLUE and RICE) are also indicated 
(from |67l). 



minating parts of the universe no other particles can. The very fact that CRs exist 
with energies greater than 10^*^ eV is the greatest motivation for UHE neutrino as- 
tronomy, since there must be a neutrino counterpart to such emission. One must 
however convince oneself that acoustic neutrino detection can compete with ex- 
isting methods and support the detection of neutrinos at the highest energies. 
Table IZT] shows a comparison of signal attenuation lengths for different methods 
of detection. Acoustic detection methods can offer the largest instrumented and 
effective volumes since the signal propagation lengths are orders of magnitude 
greater than for Cerenkov detectors. 
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Table 2.1: Attenuation lengths for different techniques 

2,6.1 "Back of the envelope'' comparison of effective volumes 
for an optical neutrino telescope and an acoustic telescope 
with 1000 sensors 

The maximal effective volume V^a of ANTARES to v^^ + iT^ is - 1 km^ |[68l. ANTARES 
comprises approximately 1000 sensors, instrumented in a volume of 0.03 km'^. If 
one assumes a regular spacing between sensors then the volume occupied per 
sensor is approximately 0.03/1000 = 3 x lO'^km^/ sensor. The simplest arrange- 
ment of these volumes is a simple cubic lattice as illustrated in Figure I2.14[ 




Figure 2.14: The instrumented volume is divided equally amongst each sen- 
sor. In a simple cubic lattice each sensor sits at the centre of a sphere of radius 
r. 

This gives a separation between sensors of: 

-Tir = 3 X 10"^ 
r = 0.0192 km 
2r = 40 m 

which corresponds to 80% of the ~ 50 m attenuation length of 450 nm Cerenkov 
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light. 

The most conservative estimate for Vgff of an acoustic array is to assume that 
Vgff = ^inst where Finst is the volume instrumented by the detector. If one popu- 
lates a volume equivalent to 1000 acoustic sensors separated by 80% of an again 
conservative 10 km attenuation length for sound: 



3 1000 
^eff = 1000 X ^7r(8/2)3 
Veff = 85333km3 



Remember, we have already been conservative in our approach, and have 
shown that if Vgff = V^inst for ari acoustic sensor array, it can have an effective 
volume ^ 85000 times greater than an optical array, with the same number of 
sensors. It has, however, been suggested that an acoustic array can experience 
Vgff up to one hundred times greater than \4nst II69IL which would then give a 
85 X 10^ times larger effective volume for an acoustic array over an optical. 



caveat The Optical Modules in ANTARES are not spaced on a regular cubic lat- 
tice, the exact geometry is complex and is dependent on reconstruction re- 
quirements, but since the reconstruction requirements of our hypothetical 
acoustic array are unknown, we give both detectors an equivalent geome- 
try and an equivalent sensor separation as a function of signal attenuation. 
Similarly positioning acoustic sensors with a separation of 10 km, just be- 
cause of the long attenuation lengths of sound, neglects any requirement 
for reconstruction. 



In conclusion, an acoustic array can potentially have an effective volume of the 
order 85 x 10"^ ^ 85 x 10^ times greater than an optical array with an equivalent 
number of sensor elements, neglecting requirements for reconstruction. In the 
following chapters, for historical reasons, typically one thousand acoustic sen- 
sors are considered in a volume of one cubic kilometre, a considerably denser 
population than argued here. However, as the discussion progresses a sparser 
optimum sensor density, taking into consideration reconstruction requirements 
will be considered. 
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2.7 Summary 

Three methods of detecting cosmological neutrinos have been discussed. Each 
one is at a different level of maturity yet all are proving fruitful. Figure 12.151 in- 
corporates the present and predicted sensitivities from Figures |23] and IZTl into an 
illustration of the energy ranges for each method of detection. The expected en- 
ergy range of acoustic detection is also included. The case for the development of 
acoustic detection lies firstly in its sensitivity to the highest energy neutrinos, thus 
completing the study of the entire neutrino energy spectrum, and secondly, the 
ability of the acoustic technique to deliver detectors with vast effective volumes. 
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Figure 2.15: The range of neutrino energies detectable by each technique 
illustrated in terms of the existing and projected sensitivities of several neu- 
trino detection experiments. 
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Chapter 3 

Simulating Neutrino Interactions 



3.1 Introduction 

This chapter deals with the simulation of neutrinos colliding into the Earth at 
energies of several Joules. It describes the first links in the simulation chain that 
ultimately predicts the sensitivity of a hypothetical array of underwater acoustic 
sensors to a flux of ultra high energy neutrinos. The first section discusses the 
underlying event - the neutrino deep inelastic scatter, the second section describes 
simulation of the induced particle shower and the final section describes how the 
acoustic signal is generated as a result of the thermal energy deposition of the 
hadronic cascade. 



3.2 Neutrino Event Generation 

The kinematics of neutrino deep inelastic scattering were introduced in the open- 
ing paragraphs of Chapter |2l In order to determine the contribution of the neu- 
trino energy to the induced hadronic cascade it is necessary to compute the Bjorken- 
y dimensionless scaling variable. Physically y describes the fraction of the inci- 
dent neutrino energy that is carried away by the hadronic system. 

The All Neutrino Interaction Simulation (ANIS) II70II has been developed for 
analysis of data from the AMANDA p9|| neutrino telescope. It is a fully object 
oriented C++ toolkit, utilising the CLHEP vector class and HepMC event records 
|7T|. The program generates neutrinos of any flavour according to a specified 



flux and propagates them through the Earth eventually forcing them to interact 
in a specified volume should they not be attenuated en route. The propagation of 
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low energy (< lO'^GeV) neutrinos through the Earth's crust and the attenuation 
of neutrinos at high (> 10^ GeV) energies is plotted in Figure |3.1[ The range of 
neutrinos, through the Earth's crust, at low and high energy is plotted in Figure 

in 




Figure 3.1: Low energy neutrinos (Left) propagate freely from the surface of 
the Earth (+ markers), through the crust to the fiducial volume (• marker) 
where they are forced to interact. High energy neutrinos (Right) do not travel 
far from their point of origin (+ markers) before they are attenuated by the 
crust (x markers). This demonstrates the increase in the neutrino nucleon 
cross section with neutrino energy. 

It is assumed that all neutrinos interact in a fiducial volume, a cylindrical 
"can" surrounding the volume instrumented by hydrophones. The energy spec- 
trum is assumed to be flat in log£, hence dN / dE oc E~^. Only neutrinos origi- 
nating from a positive hemispherical shell (with its origin at the centre of the in- 
strumented volume) are considered. Hence the ANIS program need only propa- 
gate neutrinos through the detection medium (water) and not through the Earth's 
crust. 

The mean fraction of the neutrino energy imparted to the hadronic cascade is 
approximately 25%, however, the actual value can vary from — 100%. This is 
shown in Figure 13.31 and is assumed to be independent of neutrino energy. Ini- 
tially in this simulation all hadronic cascades are assumed to contain 25% of the 
energy of the incident neutrino, since vertex reconstruction relies on the timing 
and not the magnitude of a pressure pulse. The effect of a non constant y-value 
will be discussed in Section |5^ later when describing attempts to reconstruct the 
energy of a neutrino from a given event. There is one remaining issue to be ad- 
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Figure 3.2: Range of neutrinos in ANIS at 1 < E,, < 10^ GeV (dotted) and 
10^ < El, < 10^^ GeV (solid). The low energy neutrinos can penetrate up to 
the full diameter of the Earth ~ 12000 km before they are forced to interact, 
whereas high energy neutrinos are attenuated by about 2000 km of rock. 




50 



(a) Energy dependence 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Bjorken-y 

(b) Bjorken-y probability density for all ener- 
gies 



Figure 3.3: Bjorken-y probability density for three flavours of neutrino in 
the energy range 10^ — > 10^^ GeV, interacting in water. As indicated by plot 



3.3(a) the shape of the probability density is independent of the neutrino en- 



ergy. Plot 3.3(b) is the three flavour y-distribution averaged over all energies 
between 10^ 10^'^ GeV (i.e. a projection of plot 3.3(a) along the energy 
axis). 
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dressed at this stage of the simulation, that of event multiplicity. In the section 
that follows it is assumed that the hadronic cascade is initiated by a single excited 
hadron that, following the neutrino DIS, carries 25% of the energy of the incident 
neutrino. More accurately there will be some small number of excited partons l] 
between which this percentage of the neutrino energy is shared; however, ANIS 
restricts itself to calculation of the energy of the excited hadronic final state and 
not the parton content. 

A comparison of Bjorken-y distributions from ANIS and two other event gen- 
erators, namely PYTHIA and a bespoke program developed for the ACoRNE 
collaboration, (see ||72|) has been made. The results in Figure l334l indicate that the 
produced distributions are consistent. 




to I I 

0.1 0,2 0.3 0.4 0.5 0.6 0.7 0,8 0.9 1 

y 

Figure 3.4: Bjorken-y distributions at 10^ GeV from different event genera- 
tors. 



3.3 Simulating Neutrino Induced Particle Cascades 

3.3.1 Physics processes 

Given an expression for the energy of the hadronic system one can proceed to 
simulate the development of the subsequent particle cascade within a given ma- 
terial. A suitable particle physics toolkit for simulating the passage of particles 

^i.e. some admixture of valence quarks, sea quarks and gluons 
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through matter is GEANT4B |73| . There are three categories of physics model at 
the heart of this program: those driven by theory; parameterised models which 
combine theory and data; and empirical physics models driven purely by data. 
The user is required to choose which physics models and which cross-section 
data are used in a given energy range. Two models can overlap in energy so long 
as one model does not fully occupy the same total energy range as another. A set 
of ready made high energy calorimetry "physics lists" are distributed with the 
Geant4 source code. Conceptually the formula for construction of a physics list 
is as follows: 

physics list = particles + processes + cross-sections (3.1) 

Unless a particle is assigned a process it will do nothing in the simulation. If 
a particle is assigned more than one process then the processes compete. One 
process may invoke many models and each model has a default cross-section. It 
is first determined when and where a process should occur, which depends on 
interaction length and cross-section; secondly the final state is generated, which 
is dependent on the model invoked. 

3.3.2 Particle production thresholds 

A particle propagates through the simulated detector losing energy via the pro- 
duction of secondary particles. There must be some energy below which a par- 
ticle no longer produces secondaries otherwise the program will suffer from in- 
frared divergence. Hence the user must impose an energy threshold cut. How- 
ever, such a cut may result in a poor estimation of the stopping location and 
energy deposition of the particle, so, the cut is made on the particle's range in- 
stead. The range cut is the same for all materials but the corresponding energy 
threshold is material dependent. Once a particle reaches the energy below which 
no further secondaries are produced it is tracked to zero energy through a con- 
tinuous energy loss mechanism. 

In order to generate the thermal energy density E(r') required by the pres- 
sure field integral in Equation 12.71 one initiates a shower simulation with a single 
proton at an energy corresponding to 25% of the neutrino energy. The energy de- 
posited by each successive interaction in the simulation is recorded in a RQQT[|74|| 

^results from version 8.0 are presented throughout 
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ntuple that can be analysed offline. As the shower ages its composition tends to- 
ward electrons and gamma rays through pion decay; the range cuts chosen for 
water correspond to a 3 MeV energy threshold for tracking of e^ and 7. Below 
this energy no further secondaries are produced (a gamma for instance will no 
longer pair-produce) and the particle is forced to deposit the remainder of its en- 
ergy continuously. This dramatically reduces CPU time and prevents the ntuple 
file size from diverging. The plots in Figure |331 illustrate the shape of the energy 





(b) Longitudinal Energy Deposition 



Figure 3.5: Average radial (Figure 3.5(a)| ) and longitudinal (Figure 3.5(a) I en- 
ergy profiles for one hundred 100 TeV Geant4 proton showers with (solid) 
and without (dotted) range cuts. 
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deposition with or without cuts applied. Whilst the effect on the shower shapes is 
negligible when compared to individual shower fluctuations (on top of which the 
average distributions are superimposed) the CPU time per event decreases from 
O 20 hours to O 20 minutes at 100 TeV. Furthermore the low energy extremes of 
the shower (where the agreement is less good) have little effect on the resulting 
acoustic pulse, the shower core being the dominant contributor. 

Four theory driven high energy hadronic calorimetry "use-by-case" physics 
lists are provided with the GEANT4 distribution. They are composed of the Quark 
Gluon String (QGS) model or the FRITIOF (Lund string dynamics) model [75l 
combined with either Pre-Equilibrium or Chiral Invariant Phase-space (CHIPS) 
decay modes. The range of validity for each model extends from lOGeV to 
100 TeV. The effect of each physics list on the radial and longitudinal shower 
shapes is shown in Figures 13.61 and 13.71 respectively. The CHIPS decay mode 
tends to produce a broader, older shower, with a central energy "hole" in com- 
parison to the pre-equilibrium decays, whilst there is little to distinguish between 
the QGS and FRITIOF hadron interaction schemes. References for these models 
can be found, for example in HZH, IZZll, ESl and EUl, as cited in the Geant4 doc- 
umentation. The Pre-Equilibrium model is favoured over CHIPS since it is more 
consistent with the distributions seen in UHECR air showers and it is the model 
preferred by other groups undertaking acoustic studies. Comparisons of the lon- 
gitudinal and radial shower shapes produced by GEANT4 and a modified version 
of the CORSIKA air shower program are plotted in Appendix |Bl 

3.4 Formation of the Acoustic Signal 

The acoustic signal resulting from a hadronic cascade is computed by integrating 
the energy contained in the ROOTntuple as produced by the GEANT4 simula- 
tion according to Equation 12.71 A typical pulse resulting from a 100 TeV proton 
induced shower is plotted in Figure 13. 8[ 

For the purpose of simulating the performance of large scale hydrophone ar- 
rays it is necessary to determine a relationship between the energy of a neutrino 
and the resulting pressure amplitude at a given location. The shape of the cascade 
energy density becomes more uniform with increasing energy. The integrated en- 
ergy of the cascade scales linearly with the energy of the neutrino; hence, the mag- 
nitude of the resultant pressure signal is assumed to scale with the energy of the 
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Figure 3.6: Average radial energy deposited per 20 g cm^-^ vertical slice per 
unit radial distance for one hundred proton showers at 100 TeV, using each 
of the high energy hadronic calorimetry physics lists. Depth refers to the 
distance along the shower axis. 
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Figure 3.7: Average longitudinal energy deposition of one hundred 100 TeV 
proton showers, at various depths along the longitudinal axis of the shower, 
using each of the high energy hadronic calorimetry physics lists. 
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Figure 3.8: The pressure pulse for a 100 TeV proton induced Ge ant4 shower 
in seawater, at a distance of 1km from the cascade. The angle with the plane 
transverse to the shower direction is ° . 



neutrino. Consequently an analytical parameterisation is sought to relate peak 
pressure amplitude to neutrino energy. The details of how this parameterised 
pressure is attenuated as it propagates from the cascade to a given hydrophone 
is discussed in the next chapter. The peak pressures at 1 km from various proton 
induced GEANT4 cascades at energies in the range 10"^ 10^ GeV are plotted 
as a function of the proton energy in Figure |3.9[ A linear function is fitted to 
the values allowing for a direct determination of the peak pressure at 1 km from 
an arbitrary shower for a neutrino of a given energy. The gradient indicates a 
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mean pulse peak pressure of 2.12 x 10~^^ Pa per GeV of thermoacoustic energy. 
The subsequent pulse at each hydrophone location is suitably scaled and then the 
effects of attenuation are applied. 
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Figure 3.9: Plot of mean peak pulse pressure against energy of the hadronic 
shower demonstrating the linear relationship between the two. Error bars 
are of magnitude Idmean- The residuals plot shows the deviation of the data 
points from the fitted function, there is an inherent instability at higher en- 
ergy due to the limited accuracies of the models used (see Appendix|B]|. 



3.5 Summary 

From the simulation of UHE neutrino interactions we have learnt that approxi- 
mately 25% of the energy of a neutrino is carried away by the hadronic compo- 
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nent of neutrino-nucleon deep inelastic scatter. It is therefore assumed that the 
particle cascade, and subsequent thermoacoustic energy deposition, induced by 
a neutrino with energy Ey, is equivalent to that produced by a proton at £i//4. 
Furthermore, from the simulation of the evolution of TeV proton showers in the 
sea, we have learnt that the peak pressure of an acoustic pulse resulting from 
such a shower at a given distance, scales linearly with the energy contained in 
the shower. These findings equip us with a parameterisation allowing us to esti- 
mate the peak pressure p of a pulse at 1 km from a particle cascade initiated by a 
neutrino of energy Ey given by: 

p = 0.25 xkyXEy (3.2) 

where ky is the conversion constant and is equal to 21.15 x 10~^^PaeV~^. The 
parameterisation in Equation 13.21 will be used in the remainder of the simulation 
chain discussed in the chapters that follow. 
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Simulating the Propagation of the 
Acoustic Signal 

4.1 Introduction 

Thus far it has been reported how the mechanisms which give rise to the distinc- 
tive bipolar pressure pulse, expected following an UHE neutrino interaction in 
the sea, are computed. This Chapter subsequently introduces parameterisations 
of signal attenuation resulting from its propagation through seawater. Incorpo- 
rated into this attenuation are the effects of refraction, resulting from variation in 
sound velocity in the sea, and the effects of ambient noise and filtering. 

4.2 Attenuation of the Acoustic Signal 

A hydroacoustic pressure wave is attenuated as it propagates away from its point 
of origin. There are three contributions to this attenuation: a geometric factor 
whereby the pulse diminishes the farther it travels, a diffractive attenuation that 
is a consequence of the coherent nature of the acoustic emission, and the attenu- 
ation of the pulse due to the bulk mechanical properties of the medium through 
which it propagates. Each of these effects is discussed below. 

4.2.1 Geometric attenuation 

As described in Section 12.5.21 the development of the particle cascade occurs in- 
stantaneously on the thermal diffusion timescale. The result of this is coherent 
emission along the shower axis. In the far-field (r ^ I where r is the distance 
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from the interaction and I is the cascade length) the geometric attenuation of such 
emission goes as the reciprocal of the distance travelled 1159 II . 

4.2.2 Attenuation by the medium 

The second contribution to the attenuation of the acoustic radiation, produced by 
a neutrino induced particle shower, is governed by the bulk mechanical proper- 
ties of the medium through which it traverses. Water is a low-pass filter, high 
frequency components being readily attenuated, this can be seen in the sound 
attenuation co-efficient plot from Lehtinen et al [62| shown in Figure 14.11 For the 
purposes of the large scale detector simulation the attenuation was parameterised 
by a fourth order polynomial and recomputed for a Mediterranearj^ water tem- 
perature of 25 °C. The effect of this attenuation on a bipolar signal is shown in 
Figure |42l 



10^ 




Frequency (Hz) 



Figure 4.1: The sound attenuation coefficient for seawater (from ||62|). Also 
shown is the absorption coefficient for distilled water. 

The primary reason for the difference in absorption coefficient between sea- 
water and distilled water is the excitation of magnesium sulfate and boric acid 
molecules that are present in seawater. Such molecules can have multiple sta- 



The Tyrrenian Sea located north of Sicily was used in the ACoRNE funding proposal llSOl 
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ble states and transitions between states can be induced by a change of pressure, 
effectively converting hydroacoustic pressure into heat. 
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Figure 4.2: Variation of pulse shape with respect to distance, due to seawater 
absorption, for a signal normalised to 1 Pa at 100 m from the particle cascade. 



4,2.3 Angular spread 

As one moves out of the plane perpendicular to the shower axis, that is, the plane 
of the pancake, the acoustic radiation, and hence pulse pressure, is attenuated 
due to diffractive effects. This frequency dependent process both diminishes the 
pulse height and broadens the pulse width. For each frequency component of 
the attenuated pulse, there is a corresponding sin(x)/x diffraction pattern, in 
analogy to the diffraction of light through a narrow slit. The length of the shower 
corresponds to the size of the aperture. The pulse at zero degrees out of the 
pancake plane is computed by taking the Inverse Fourier Transform (IFT) at zero 
angle. The pulse at a given angle out of the pancake plane is thus computed by 
taking the IFT of the respective sin(x) / x function for that angle IISTl . This effect is 
plotted in Figure |43l The pulse amplitude is attenuated to nearly one thousandth 
of its original magnitude at 5 ° out of the pancake plane for a 10 m long shower. 
The angular distribution plotted in Figure 14.31 is parameterised and then used to 

■'httpV/www.mathworks.com/ Accessed 02/05/2006 
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Figure 4.3: Variation of far-field pulse shapes (Left) and signal amplitudes 
(Right) as a function of angle out of the plane of the pancake (from HH). 
The step in angle of 0.26 ° derives from the digitisation of the sample into 
1024 bins, which was chosen because it facilitates a faster IFT algorithm in 

matlabB 



scale the peak pressure recorded on a given receiver as a function of its angle out 
of the plane perpendicular to the shower. 

The three components of attenuation described above are applied in succes- 
sion to the unattenuated peak pressure amplitude as calculated via the linear 
energy to pressure relation formulated in Section I3.4[ Given a formalism for 
calculating the attenuation of an acoustic pulse, as a function of the distance it 
propagates through the sea, a method for determining the distance travelled is 
desired. 



4.3 Refraction 

The velocity of sound in seawater is dependent on temperature, pressure, and 
salinity. A simple means to describe the stratification of the sea into horizon- 
tal layers of constant sound gradient is illustrated in Figure 14.41 There are three 
zones of sound propagation to consider: near the surface there exists an isother- 
mal layer which is well mixed by surface waves and wind, within the isotherm 
the sound velocity increases with depth due to a slight increase in pressure; be- 
low the isotherm is a region known as the thermocline through which the sound 
velocity decreases with depth due to the corresponding decrease in temperature; 
below the base of the thermocline, in the deep ocean zone, the sound velocity gra- 
dient inverts and sound velocity increases with depth due to the corresponding 
increase in pressure. Sound rays tend toward the region of lowest sound veloc- 
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Figure 4.4: Three layers representing regions of constant sound gradient in 
the sea: the isotherm, the thermocline and the deep. 

ity. Where there is an inversion in the sound velocity gradient sound rays are 
refracted away from each other, leading to the formation of "shadow zones", re- 
gions from where acoustic signals may be inaudible. When building an acoustic 
neutrino telescope it is desirable to be far away from such environments since 
this will contribute to a loss in effective volume. A full derivation of the ray trac- 




z 



Figure 4.5: Geometry of a single curved ray in a medium with varying 
sound velocity. 

ing equations in the general case, where the sound velocity gradient is not linear, 
and for the case of a linear sound velocity gradient can be found in Appendix ICl 
For the geometry illustrated in Figure 1431 with a constant sound velocity gradient 
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g = dc{z) I dz, it can be shown that: 



Xi — Xq 



Co / sin ^0 — sin di 



) 



(4.1) 



g V COS^o 



g COS ^0 



(4.2) 
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g COS ^0 



(cos 01 — COS 00 ) 



(4.3) 



h - k 




cos 0o(l + stn0i) 



} 



(4.4) 



cos 01 (1 + stn0o) 



where Co = c(z = 0) is the sound velocity at the origin depth, is the time at the 
point of origin and ti is the time at point (xi,zi). Angles 6q and 0i are the angles 
the ray makes with the horizontal at points {xq,zo) and (xi,Zi) respectively; Sq 
and si define the arc length along the ray. For simplicity the problem has been 
reduced into two dimensions defining the xz-plane (since refraction only occurs 
in the vertical (z) direction one can arbitrarily align rays along one of the hor- 
izontal axes). The ray tracing algorithm iterates in steps of x for a given path 
between source (cascade) and receiver (hydrophone). Equation 14.21 can be solved 
for the deflection angle 0i which then forms the input to Equation l4.3l along with 
00 which is known. 

In order to compute the trajectory of the real refracted ray R that intersects 
source and receiver at {xq, zq) and (xi, Zi) respectively one must first compute the 
angle of deflection 0^ of the deflected ray d that results from sending a ray along 
the imaginary linear path i. This is because only the sound velocity as a function 
of depth, the coordinates {xq,zq) and (xi,zi); and the angle 0, are known a priori. 

A plot of the ray trace as computed in the large scale detector simulation as 
described in Chapter |5] is shown in Figure 14.71 In the case of a constant sound 
velocity gradient: 



where c is the sound speed at depth z and Cq is the sound speed at z = 0. The 
values for cq and g are obtained from a sound velocity profile (SVP) for a partic- 
ular location in the sea. The SVP at the site of the ANTARES neutrino telescope 



g = constant 



c = 



co + gz 
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Figure 4.6: The difference between the angle of deflection 6^ of the deflected 
ray d and the start angle 0/ of the undeflected imaginary ray i is A9 = 9d—0j. 
The start angle for the real refracted ray R is therefore given by 0j^=0,— A0. 



is shown in Figure 14.81 superimposed on a ray trace for a source emitting from a 
depth of 1500 m at various angles to the horizontal. The effect of the inversion 
in the sound velocity gradient at the base of the thermocline is included, hence 
rays are refracted back toward the sea bed as they approach the surface. In the 
hydrophone array simulation only the deep sea zone, with approximated linear 
SVP, is considered, such that there is no inversion in the sound velocity gradient 
and rays are terminated at the surface without reflection. Similarly rays that hit 
the sea floor are terminated without reflection. This represents a more conser- 
vative approach than one where boundary reflections are considered, since the 
effective volume will increase if rays can reach a given array via surface or sea 
floor reflections. Detailed information of the boundary conditions is required to 
effectively model ray reflection, hence it is presently neglected. In the simulation 
the hydrophone array sits at —2 km in the deep ocean zone, far away from the 
thermocline in a region where the sound velocity gradient is linear. The effect 
of a linear sound velocity profile is to send acoustic rays along circular trajec- 
tories H82]| of radius R = co/g = 1506/0.0156 ^ 93 km. There exists a natural 
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Figure 4.7: A ray trace from the large scale detector simulation. The solid 
curved ray and the dashed curved ray correspond to the rays R and d re- 
spectively in Figure l46l 

upper limit on source distance, due to refraction, in a sea of finite depth and con- 
stant sound velocity gradient. The concept is illustrated in Figure |4.9t sources 
greater than 48 km beyond an array located at mid depth in a 4000 m deep sea 
cannot be detected. This is the absolute maximum distance from which a source 
can be heard due to the geometrical constraints of refraction. With respect to the 
performance of an acoustic array, by demanding several signals in coincidence, 
which are themselves emitted subject to angular constraints, it is unlikely a neu- 
trino induced particle cascade could be successfully reconstructed at a distance 
of this magnitude. The exact limit on the distance from which cascades can be 
reconstructed is of importance when calculating the effective volume of a given 
hydrophone array and will be discussed further in Section |53I 

Successful computation of curved ray trajectories thus allows the calculation 
of hit times and ray path lengths from Equations 14.21 and 14.41 respectively. The 
effect of attenuation as discussed in Sections 14.2.11 to 14.2.21 is thus applied to the 
refracted times and ray paths. 

4.4 Ambient Noise, Filtering and Thresholds 

Estimation of the sensitivity of a hydrophone array relies on reconstruction of 
the neutrino interaction vertex from the acoustic signals received. Since it is the 
timing of the signals that is of importance in the detector simulation, a parame- 
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Figure 4.8: Ray trace for a source at 1500 m depth. Refracted rays (solid) 
and unrefracted rays (dashed) are shown. The SVP for the ANTARES site is 
also shown, with sound speed across the upper x axis. The ANTARES SVP 
is essentially linear. Neglecting the inversion in the sound speed gradient 
equates to extrapolating the linear section directly to the upper x-axis (This 
yields Co = 1506 ms^^ andg = (-)0.0156s-i). 



terisation of neutrino energy to peak pressure can be used. In reality event trig- 
gers will depend on knowledge of the bipolar pulse shape. The performance of 
a matched filter is coupled to the attenuation calculation introduced in Section 
14.2.31 as a means of simulating the response of a hydrophone to bipolar signals. 
The matched filter is applied to a flat, Gaussian white noise background. The 
noise level is determined from Rnudsen curves (see Appendix I A.2|) which relate 
the surface wind speeds to the level of ambient noise. The matched filter en- 
hances the signal by approximately a factor of three. The simulated hydrophones 
are modelled to have omnidirectional sensitivity; an example of matched filter 
performance is plotted in Figure I4.10[ Naturally the precise shape of a bipolar 
signal in the sea varies as a result of the attenuation described in Section l4!2i As 
such a particular matched filter will be optimal for some function of distance and 
angle from a given source. In reality several matched filters will have to be used 
to fully optimise neutrino signal identification. This ensemble of filters can be 
represented in the large scale simulation by a single entity, optimised for all fre- 
quencies, with a probability for correct identification of a neutrino induced signal 
equal to unity. 

The signal detection threshold, in the detector simulation, is set at 35mPa, 
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which corresponds to a probability of false alarm (PFA) of one false trigger in ten 
years, due to noise exclusively resulting from surface wind speed, with a five-fold 
coincidence. A threshold of 35 mPa indicates a mean noise level of 35.9 dB. The 
mean noise level from June to August in the Tyrrenian Sea for example is 35 dB 
with an average wind speed of 5.6 knots liSOi . 

4.5 Summary 

We have seen that acoustic radiation is modified by its transmission through the 
sea. This leads to a set of constraints on the simulation of a large scale acous- 
tic array. Hydrophones that are greater than 5 ° out of the plane defined by the 
acoustic pancake are neglected because of the signal attenuation resulting from 
diffraction. There is a natural, geometry dependent, limit on the distance from 
which a source can be heard by a hydrophone located in a sea of finite depth be- 
cause of the refraction of sound rays towards the surface. Finally, a threshold for 
detection of 35 mPa has been established based on the ambient noise in the sea 
and the probability that this noise could emulate a neutrino like acoustic signal. 



66 



4.5. Summary 




''^ ^ sea bed 

dx 



Figure 4.9: Rays following circular trajectories of radius R = 93 km and 
originating from the shaded zone and beyond cannot reach the array. Hence 
there is a natural distance D = 2dx beyond which sources cannot be heard 
(so long as rays are not reflected and any inversion in the SVP is neglected, as 
is the case here). The Mediterranean has a maximal depth of about 4 km, ide- 
ally an array would sit equidistant between the base of the thermocline and 
the sea bed. In this simulation the isotherm and thermocline are neglected 
along with surface and sea-bed reflections, and the array sits at mid depth. 
dx = y/R^ - (R - dy)^, yielding a value for D = 2 x = 2 x 24 = 48 km. 
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20nnPa Peak Bipolar Pulse Injected at n =0 
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Original Data 
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Figure 4.10: The effect of a matched filter on a bipolar signal embedded in 
some Rona data (from ||83| which contains a full description of the matched 
filtering technique). The effect of such a filter on a signal embedded in a flat, 
Gaussian, white noise background is to enhance signal to noise by approxi- 
mately a factor of three; this is the case in the large scale detector simulation. 
The lower plot is the raw input with the bipolar signal emboldened because 
it is difficult to distinguish by eye from the noise. Offset above is the filtered 
output, where the bipolar signal is clear relative to the background noise. 
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Large Scale Detector Simulation 



5.1 Introduction 



All the ingredients are now available for the recipe that computes the sensitivity 
of a hypothetical array of hydrophones to a flux of UHE neutrinos. In summary 
the following steps have been described: 

• Generation of the mean fraction of the incident neutrino energy taken away 
by the hadronic particle cascade. 

• Simulation of the evolution of a hadronic cascade through seawater, yield- 
ing the thermal energy density deposited therein. 

• Integration of the cascade energy and subsequent production of the bipolar 
acoustic signal. 

• Parameterisation of the pulse formation process for fast calculation of peak 



acoustic pressure using the far-field approximatiori^. 

• Modelling of the attenuation of acoustic radiation as it permeates the sur- 
rounding seawater thus diminishing the peak pressure. 

• Calculation of the trajectories and signal times for sound rays undergoing 
refraction due to a gradient in the sound velocity as a function of ray depth. 

The following discussion explains the procedures undertaken in a bespoke simu- 
lation for the purpose of estimating hypothetical hydrophone array performances. 

near-field pulse will obviously contain much higher frequency components, but the as- 
sumption that the matched filter is optimal across all frequencies of interest takes this into con- 
sideration 
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5.2 The Sensitivity Estimation Procedure 

The steps taken by the large scale detector simulation are as follows: 

Step 1. An ensemble of downward going neutrinos is generated with an energy 
spectrum linear in log-|^Q Ey. 

Step 2. A hydrophone array is "constructed" in a deep sea environment from a 
file containing the coordinates of typically one thousand hydrophones dis- 
tributed at random in a volume of one cubic kilometre. 

Step 3. The program then samples the neutrino spectrum and calculates the de- 
tector response on an event by event basis: 

1. Each neutrino is forced to interact in a fiducial volume, called the "can", 
that surrounds the volume defined by the hydrophone locations. 

2. The peak pressure as a function of neutrino energy and attenuation 
at each hydrophone is calculated for both refracted and unrefracted 
sound rays assuming an omnidirectional sensitivity. 

3. The pressure amplitude and arrival time for all those hydrophones that 
register a signal above threshold is recorded, taking into consideration 
the performance of a matched filter. 

Step 4. For those events that register hits above threshold reconstruction of the 
interaction vertex and the neutrino trajectory is attempted. 

Step 5. If the vertex reconstruction algorithm returns an unphysical vertex loca- 
tion the event is discarded, otherwise the event is classified as detected. 

Step 6. The sensitivity of the hydrophone array under examination is then calcu- 
lated from the number of successfully reconstructed events. This is trans- 
lated into a limit on the neutrino flux assuming that the hydrophone array 
detected no events for a given period of observation. 

An example of the source distributions from which the event data are sampled 
can be found in Appendix O Since we are operating close to threshold, a hard 
cut can be placed on the angle of a given receiver with respect to the plane of 
the pancake such that above 5 ° or so the receivers are neglected since the angular 
attenuation has diminished the signal by several orders of magnitude (see Section 

MM- 
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5.3 Neutrino Spectra and Geometries Used 

As with the neutrino event generation program discussed in Section |3^ the neu- 
trino energy spectrum is arbitrarily modelled to be flat in log E since the ensu- 
ing sensitivity calculation will incorporate a model independent neutrino flux. 
Therefore: 

constant 4^ -— cx E (5.1) 



rf(log Ey) dE 
The neutrinos are isotropic over a hemispherical shell that is centred at the origin 
of the instrumented volume. All neutrinos propagate freely from the shell to a 
point in the cylindrical, fiducial volume, known as the "can" that envelops the 
volume defined by the hydrophone locations. The distribution of generated neu- 
trino events is plotted in Figure 15.11 The origin of the coordinate system that is 




Figure 5.1: The spatial distribution of neutrino source locations isotropic 
over a positive hemispherical shell centred on the instrumented volume. 
Surrounding the instrumented volume is the cylindrical can. The radius of 
the shell is somewhat arbitrary as long as the shell itself encompasses the 
can. In this plot it is 100 km. 

employed in this simulation is defined to be at the centre of the instrumented vol- 
ume at a depth of 2000 m in a sea of depth 4000 m. A conventional right-handed 
set of axes is fixed thus, with the z-axis running vertically and the xy-plane defin- 
ing the horizontal. This coordinate system is illustrated in Figure 15.21 Typically 
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1000 hydrophones are placed, at random, in a volume of one kilometre cubed, 
centred at the origin. 
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Figure 5.2: Coordinate system employed in the hydrophone array sensitiv- 
ity estimation program. 



5.4 Event Reconstruction 

Signal arrival times are used to resolve interaction vertices. The formula for ver- 
tex reconstruction of sound rays that have undergone refraction is different from 
the formula that is used for rays that have linear trajectories. 



5.4.1 Vertex location with linear rays 

In a homogeneous medium the velocity of sound c is constant. Point source lo- 
cation from differences in signal arrival times is well documented and has been 
implemented by the RICE (Si, SAUND |[6Zl, Erlangen BH and ACoRNE groups. 
The standard formalism is described below. 

A receiver at location r,, hears a signal at time ti, from a source at location v, 
originating at time t^, according to the following equation 

\v-r,\^ = ^{k-tQf (5.2) 
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Expanding the bracketed terms and solving for v gives: 



(,2 .2^ .2r... ..^2 



r-)-c%U-tQY-{tj-tQy 



c^iu - toy 



2V ■ [fi - Tj) 



(rf-rp-c'ltf-tj-ltoit.-tj)] = 2v.{n-fj 



c^[tj-tj]+2c%{t,-tj) = 2v-{r,-fj 



which can be expressed in the form of a matrix equation: 

R + to^ = Mv 
where, in the limiting case of four receivers: 



R 

M 



r] - r] - ^{t] - tj) = 1,2; 1,3; 1,4) 
2c\t,-tj) (z,; = l,2;l,3;l,4) 

^ dxi2 dy\2 d.z\2 ^ 

dXi3 dl/13 rfzi3 

dx\^ dy\^ dz\^ 



(5.3) 



(Here, the receiver with index 1 is taken as the reference receiver, dx, dy, dz are 
the difference in x, y and z coordinates respectively, between this receiver and the 
other receivers that detect a signal)B Therefore: 



V 



One defines: 

y = M~^K (5.4) 
Q = M-^T^ (5.5) 

then: 

y = y + io Q (5.6) 

^it now becomes apparent why a random spatial distribution of hydrophones is employed, 
since a regular sensor array could lead to singular (| M | = 0) or near singular matrices, for example 
if the hits were confined to a plane aligned with either of the x,y or z axes; such a matrix would 
fail the inversion process 
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The minimum number of receivers required for vertex reconstruction to be pos- 
sible is four. If there are n receivers the system becomes over constrained and 
matrix M becomes a 3 x (n — 1) matrix. Because M is no longer symmetric it 
must be inverted using singular value decomposition (SVD). For an (m x n) ma- 
trix A with m > n, the singular value decomposition is given by: 

A = ULV^ 

where U and V are orthogonal m x m and n x n matrices containing the column 
and row eigenvectors respectively. L is an m x n matrix containing the eigenval- 
ues in increasing order. The inverse of A is therefore given by: 

A"i = VL"^U^ 

Since P and Q are known, it only remains to find tQ. This is done by solving the 
propagation time equation for any receiver amongst the four: 

\v-r,\^ = c^{t,-tof = 1^ + ^oQ -nf (5.7) 

which is simply a quadratic in to: 

ch} + chl - IchitQ = p2 + f2Q2 ^r^^2tQV -Q-IV -Ti- ItoQ ■ r, (5.8) 

rearranging for to gives: 

(c^ - q2)^2 + (2Q . _ 2^ • Q - 2c%)to - (p2 + rj-2'F- r, - ch]) = (5.9) 

Both roots to Equation 15.81 are potentially viable. The correct root is sought by 
minimising the following chi-squared: 

^r eceiv ers 

X^= ^^-^rf (5.10) 

where as before t[ is the arrival time of the signal at the receiver located at r,, and 
the reconstructed time ty is given, in terms of the signal emission time t±s and 
vertex location v±, from the positive and negative solutions, by: 

, , , ^i-v±) 

tr = t±s H 

C 

which yields the correct root 99.9% of the time. Interestingly, simply invoking 
causality is not sufficient since in many cases both roots are allowed; the exact 
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Mean 9.165e-07 




X resolution (m) reconstructed x(m) 

Figure 5.3: The resolution of ttie vertex x-coordinate, in tlie case of a constant 
sound velocity and no smearing on the hydrophone location, has a mean 
value consistent with zero and a standard deviatiorB^f less than 1 mm. 

number of causal and non-causal roots is dependent on the geometry of the detec- 
tor and the can. The vertex resolution of a one thousand hydrophone array, of one 
cubic kilometre instrumented volume in the case of a constant sound velocity, is 
plotted in Figure l53l Now, given the interaction vertex, and a set of receiver sig- 
nals confined to a plane perpendicular to the direction of the neutrino, all one has 
to do to resolve the cosmic source of the neutrino is point back through the vertex 
with a vector normal to the pancake plane. The neutrino pointing vector has been 
reconstructed in two ways: firstly by taking the normal to the plane defined by 
the vertex and the two most energetic receivers in the pancake, as illustrated in 
Figure l5!4l secondly by taking the cross product of two vectors orthogonal to each 
other and the cascade direction, as plotted directly from the simulation in Figure 
15.41 The procedure, subsequently called the "Top2" algorithm, for determining 
these vectors is as follows: 

1. calculate the centre of mass of the system of hit receivers - this lies in the 
plane of the pancake; then, translate the coordinate system such that the 
origin is at the centre of mass. 

2. Define a vector vi that points from the centre of mass to the vertex. 

3. Rotate the coordinate system such that Vi is aligned along (say) the x-axis, 
then collapse the system into 2D by discarding the x-coordinates. 

4. Perform a least squares fit to the hit receivers in yz-plane, because the centre 
of gravity is at the origin, the intercept of this line will be at zero, define a 

^for historical reasons the standard deviation in ROOTis actually called "RMS", an explanation 
is given at http://root.cern.ch/root/html/THl.html#THl:GetRMS Accessed 25/11/2006 
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Figure 5.4: The vector normal to the plane defined by a pair of hydrophones 
lying in the acoustic pancake and the reconstructed vertex points back along 
the axis of the hadronic cascade. Hence this vector is coUinear with the tra- 
jectory of the incident neutrino. The pointing resolution is defined by the 
angle a. between the pointing vector and the incoming neutrino trajectory. 

second vector V2 in the direction of this line. 

5. Rotate the system back - vectors V\ and V2 are orthogonal to each other, the 
neutrino pointing vector is simply Vi x V2- 




Figure 5.5: Output from the simulation displaying vectors V\ and V2 and the 
resultant pointing vector as described in the text. 

The pointing resolution is defined as the angle oc between the real neutrino 
trajectory and the reconstructed pointing vector. A comparison of the resolution 
from both methods discussed above is plotted in Figure |5^ a mean pointing res- 
olution of 6.9 ° is obtained from the Top2 algorithm and 12.0 ° from the collapsed 
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mean algorithm, in the absence of refraction, for an array of 1000 hydrophones 
distributed at random in a volume of one cubic kilometre. 



o 10^ 



Top2 mean = 6.9n 
: CO I tapsed mean = 1 2° 




50 60 70 
resolution (°) 



Figure 5.6: The pointing resolution of both reconstruction techniques for an 
array of one thousand hydrophones distributed randomly in a volume of one 
kilometre cubed and a can radius of 10 km. The resolution is defined as the 
angle between the true neutrino trajectory and the reconstructed pointing 
vector. 



5.4.2 Vertex location with curvilinear rays 

If the velocity of sound is a function of depth, c = c{z), then Equation 15.21 is no 
longer valid. Sound rays follow curved trajectories and the pancake becomes 
warped, this is illustrated in Figure 15.71 Unlike for constant sound velocity, there 
is no analytical solution to the event vertex v and the signal emission time tg. One 
now relies on iterative minimisation to interpolate towards event vertices in the 
presence of refraction. This concept has been employed elsewhere by the RICE 
||55| and SAUND |[67| collaborations to enable reconstruction of vertices in the 
data from their respective experiments. For a given set of hits on a sensor array 
one can minimise the following metric: 

^r eceiv ers 

M= £ {t,-t,f (5.11) 

where ti is the arrival time of the signal recorded by the receiver located at Fj and 
tv is the time it takes for a signal from a test location v to reach the same receiver. 
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Figure 5.7: Schematic of the energy deposition. Refraction occurs only in 
the vertical plane; there are two extreme cases: that when the shower is hor- 
izontal (Left) and that when the shower is vertically inclined (Right). In the 
former case, there is no deformation of the pancake such that the radiation 
field is bent outside the plane orthogonal to the shower axis. 

In references ||55| and |[67| minimisation begins from a location determined by 
analysis of precomputed look-up tables containing the sensor hit times for reg- 
ularly spaced locations in the fiducial volume surrounding the detector. In this 
approach, the number of test locations, and hence computations, scales with the 
size of the fiducial volume under consideration. When volumes of a hundred 
or so cubic kilometres are being dealt with the process becomes inefficient, taking 
several minutes per event (depending on can size and lattice density). Instead, an 
estimate of the vertex location is sought from which the minimisation procedure 
can begin. 

The equations for point source location in a homogeneous medium, from sig- 
nal arrival times, can be used to approximate the location of a point source in an 
inhomogeneous medium. Starting with Equation 15.31 where, in the limiting case 
of four receivers: 



Again, the receiver with index 1 is taken as the reference receiver, dx, dy, dz are 
the difference in x, y and z coordinates respectively, between this receiver and the 



R 



1,2; 1,3; 1,4) 
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other receivers that detect a signal. Additionally c = c{z), where z is the mean 
depth of the four receivers and M is always a 3 x 3 square matrix. Solving for the 
vertex v gives: 

which is once more solved via the propagation time equation: 

The mean vertex v from Nq quadruplets of receiver signals is thus calculated as: 

v = ^ E ^^ (5-12) 

There are: 

unique receiver quadruplets from N receiver signals. This number rapidly di- 
verges as N increases; for example, there are 230300 possible four-receiver com- 
binations of 50 receivers. Furthermore calculation of the unique permutations is 
in itself computationally intensive. As such a "reduced mean" with fewer ele- 
ments is desired. Three "reduced mean vertices" have been considered. The first 
is the mean vertex from each receiver and its three spatial-nearest-neighbours; 
the second is the mean vertex from each receiver and its three temporal-nearest- 
neighbours; and the third is the mean vertex from each receiver and its three 
pressural-nearest-neighbour^. In each case there are only N — 3 contributions 
to the mean, so for 50 receivers there are only 46 permutations compared to the 
230300 above. All three methods rely on Equation |5.10l to solve for the vertex, util- 
ising the sound velocity at the mean depth of a given receiver quadruplet. Now, 
neither root of the quadratic is in fact a unique solution to the propagation time 
equation (Equation |5.8|) since we have substituted a linearly variable sound veloc- 
ity in place of a constant sound velocity. As such there is an inherent instability 
in the estimation method, however, an estimate of the mean is found nonethe- 
less. Calculation of estimated vertices takes less than one second per event on a 
standard 2 GHz processor with 512 MB of RAM. 

The comparative vertex x-coordinate resolution from the three reduced mean 
vertices is plotted in Figure 15.81 indicating that the temporal-nearest-neighbours 



^i.e. if the hydrophones are sorted into a list of descending pressure, the nth phone has pres- 
sural nearest neighbours n + 1, n + 2 and n + 3 



79 



Chapter 5. Large Scale Detector Simulation 



gives the best estimate of the vertex. The x resolution of the temporal mean with 
respect to the unreduced mean is subsequently plotted in Figure 15.91 confirm- 
ing that the vertex estimate from this method is not only a computationally less 
intensive than taking the mean vertex from all receiver combinations but more 
accurate. The standard deviation on the shower time, at, for each of the reduced 
mean vertices is displayed in Table l5Jl here it can be seen why the reduced tem- 
poral mean gives the best performance, showing the smallest intrinsic error on 
the shower time as returned by Equations 15.81 and 15.101 
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Figure 5.8: The resolution of the reconstructed vertex x-coordinate. The 
solid curve represents the temporal mean with a = 1904 m; the dotted curve 
is the pressural mean with a = 2677 m; and, the dashed curve is the spatial 
mean with a = 3576 m. No smearing on the hydrophone locations was in- 
cluded and an array of 200 hydrophones distributed at random in a volume 
of one kilometre cubed was used. 10, 000 neutrinos were generated in a can 
of radius 10 km. 



Minimisation of Equation l5.11l begins from the reduced temporal mean vertex 
and is performed by MINUIT [|84|, the de facto High Energy Physics minimisa- 
tion package. A comparison of the vertex resolution of the unreduced mean, the 
reduced temporal mean and the vertex returned by MINUIT is plotted in Figure 

Islol 

Vertex minimisation for the same 200 hydrophones as for each of the esti- 
mated vertices is evaluated in Table 15. 2[ Eighty-seven percent of temporal min- 
imised vertices are found to within 1000 m of the true vertex compared to 68% of 
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X resolution (m) 

Figure 5.9: The resolution of the reconstructed vertex x-coordinate. The 
solid curve represents the temporal mean with a = 1904 m and the dashed 
curve is the unreduced mean with a = 2381 m. 



Reduced Mean 


CTtis) 


temporal 


2.68 


spatial 


2.98 


pressural 


3.04 



Table 5.1: Standard deviations of the shower times returned by Equations 
I5.8l and r5.10l for each of the reduced mean methods in a can of radius 10 km. 

estimated vertices; furthermore, 30% of minimised vertices are reconstructed to 
within 0.1 m of the true vertex compared to 1% or less of estimated vertices. Natu- 
rally different array densities will offer different reconstruction performance and 
this will be discussed further in Section l5?6l Furthermore increasing the radius of 
the fiducial volume beyond 10 km will make the effects of refraction more pro- 
nounced. 

A note on MINUIT performance 

The performance of MINUIT is not only sensitive to the starting parameters but 
the topology of the hypersurface and has not been optimised in this study. 
The minimisation procedure begins with some estimate of the true vertex as re- 
turned by the temporal-nearest-neighbour vertexing algorithm described previ- 
ously. Additionally one must specify an estimate of the error on the associated 
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X resolution (m) 



Figure 5.10: The resolution of the reconstructed vertex x-coordinate. The 
solid curve represents the temporal mean with a = 1904 m and the dashed 
curve is the minimised x resolution with a = 1259 m. 



parameters, namely the x, y and z coordinates of the estimated vertex. Further- 
more one must decide whether to constrain the values that x, y and z can take or 
whether to leave them unconstrained. In this study the vertex coordinates were 
constrained to be within the fiducial volume, i.e. the can. This was decided to be 
so because an unconstrained fit can lead to test locations well beyond the fiducial 
volume, resulting in ray traces over thousands of kilometres. This is problematic 
because the number of iterations possible in the ray trace calculation is limited by 
the amount of physical memory in the computer on which the simulation is run- 
ning. One can arbitrarily reduce the number of iterations made but this reduces 
the number of straight line segments used to approximate the circular ray path 
and thus causes it to appear more and more like a straight line therefore reduc- 
ing the accuracy of the operation. An expression for the error on the estimated 
vertex on an event by event basis should improve the minimisation procedure; 
however, for the purpose of this study the error in x and y was set at one can 
radius and the error in z was set at one can height for every event. Additionally 
the vertex reconstruction problem in the presence of refraction is non-linear and 
as such MINUIT is not guaranteed to successfully converge on the best solution. 



82 



5.5. The Effective Volume 





Temporal 


Spatial 


Pressural 


Res. < ± 


pcfiryiatpr] 


Tninimispf] 

±1. L±± LXJ.X L±l3^\^ 




TTiinimispH 

±1. Li.± Li.J.± L± 


pcfirn^afpH 


miniTTiispH 


1000 m 


yjfj /o 


87% 


43% 


81% 


64% 


84% 

/O 


100 m 


19% 


64% 


10% 


59% 


27% 


63% 


10 m 


3% 


49% 


1.8% 


48% 
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2% 
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0.1m 


0% 
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0.08% 
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1% 


29% 


0.01m 


0% 


14% 


0% 


13% 


0.7% 


13% 



Table 5.2: Minimised vertex performance. The resolution is defined as 
Xrecon — ^reai where Xyecon and Xrefl/ are the reconstructed and real vertex x- 
coordinates respectively. Minimisation beginning from the temporal mean 
vertex can be seen to provide the best performance. Here the can radius is 
10 km. 



5.5 The Effective Volume 

The effective volume V^^^ of a simulated array is defined as: 

Veff = ^ X V (5-14) 

where Ndet and Ngen are the number of events that are detected and the number of 
events that are generated respectively, and Vgen is the volume in which events are 
generated. The behaviour of Equation 15. 141 is such that by increasing the volume 
Vgen in which event vertices are generated, V^jj will increase proportionally, how- 
ever this is offset by a reduction in the number of detections, particularly at lower 
energy and in essence therefore raises the threshold. Vg^n should certainly extend 
beyond the instrumented volume and in practice one usually fixes the size of Vgen 
with the maximum range for source detection, remembering that beyond 24 km 
parts of the sea are geometrically inaccessible to sound rays subtending the array 
(see Figure W9\ in Section I4.3)) . Running the simulation with no pressure thresh- 
old and a can radius greater than 48 km determines the extent to which events 
can be detected as constrained by geometry only. The range of those events that 
produce a reconstructable vertex is plotted in Figure 15.111 This range determines 
the radius of the volume Vgen when calculating the effective volume, V^.^^. 

Given an expression for the effective volume, with a 48 km can radius it is 
possible to test the effect of different array densities. An optimum number of 
hydrophones is sought, that is, the minimum number of sensors required above 
which there is no significant improvement in the effective volume. This relation- 
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range (km) 

Figure 5.11: Range of events accessible due to the limitations imposed by re- 
fraction. The lighter histogram corresponds to the range of events detected if 
refraction is neglected and shows a continuous acceptance up to 48 km. The 
darker histogram is the range of detected events in consideration of refrac- 
tion. Beyond 24 km the acceptance reduces until no events are detected at a 
range just below 48 km. 

ship is plotted in Figure [S^H Above a density of 400 hydrophones per cubic kilo- 
metre the gain in effective volume is small, especially above 10^^ GeV. At lower 
energies sensor density is of more importance; at least 800 hydrophones per cubic 
kilometre are required for an optimal ^^e//- The matter of hydrophone density still 
requires further discussion w^ith respect to the resolving power of a given array. 

5.6 Optimal Hydrophone Density for Reconstruction 
in the Presence of Refraction 

The sensitivity of a hypothetical detector relies on its ability to reconstruct the 
vertex of a given event. However if such a detector is to be used as a "telescope" 
for the purpose of Neutrino Astronomy, then it must also provide information on 
the trajectory and energy of a neutrino. One may therefore pose the following 
question: "what is the optimal hydrophone density required to provide the best 
vertex and pointing resolution, via the reconstruction algorithms introduced thus 
far?" Typically one thousand hydrophones have been considered in a volume 
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of one kilometre cubed. From Section 15.51 a minimum number of sensors has 
been found for which the effective volume of an array is maximal. An optimum 
sensor density for vertexing and pointing accuracy is now desired. The standard 
deviation on the vertex x, y and z coordinates as a function of hydrophone density 
and the number of reconstructions within 10 cm of the actual vertex location are 
plotted in Figure 15.131 The pointing reconstruction as a function of hydrophone 
density is evaluated in Figure I5.14[ 



5.7 The Effect of Uncertainty in Hydrophone Locations 

The ANTARES |T9| neutrino telescope utilises a High Frequency Long Base Line 
(HFLB) acoustic positioning system, operating at 40 — 60 kHz over a base of 300 x 
300 m, combined with local tilt-meter-compass sensors to obtain the location of 
its sensor elements in the depths of the Mediterranean Sea. The system provides 
spatial positioning to an accuracy of 0.1 m. One may therefore assume that any 
future undersea detector array will have at least comparable performance. The ef- 
fect of smearing the hydrophone locations by a Gaussian distribution with width 
(T^ = 0.1 m on the standard deviation of the vertex x-coordinate {Ux) and the Top2 
pointing vector (c7xop2) is summarised in Table l53l 
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Figure 5.13: Resolution (Left) and accuracy (Right) of vertexing as a function 
of the number of hydrophones. Resolution in x {• marker) and y (■ marker) 
plateaus beyond a hydrophone density of a 300 per cubic kilometre whereas 
the number of reconstructions up to a given accuracy continues to improve 
up to a hydrophone density of 1500 per cubic kilometre. Because the can is 
much larger radially (48 km) than it is in height (4 km) the z-coordinate (▲ 
marker) resolution is much better than in x and y. 





Vertexing 


Pointing 


aj, (m) 


(Tx (m) 


Cx < 0.1 m 




C^Top2 < 5 ° 


0.0 


11230 


47.6% 


28.0 


8.89% 


0.1 


11800 


47.5% 


28.6 


8.60% 



Table 5.3: The effect of a 0.1 m hydrophone location uncertainty on recon- 
struction performance of a 1000 hydrophone, 1 km'^ array in a can of radius 
48 km. 

5.8 Reconstruction of the Neutrino Energy 

As previously stated, an estimate of the energy of an interacting neutrino is re- 
quired if one is to use a hydrophone array for neutrino astronomy. To be able to 
identify a source in the sky is one thing, but to understand the processes occur- 
ring therein some kind of calorimetry is necessary. A simple method was tested 
for determining the energy resolution. The following was evaluated: 

-^hydros 

X^= Ph-Pe (5.15) 

h=l 

where Nhydros is the number of hydrophones in the array, is the pressure 
recorded on each hit hydrophone and Pg is the pressure expected at each phone. 
The pressure expected at each phone is a function of the reconstructed vertex 
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Figure 5.14: Top2 (■ marker) and Collapsed (▲ marker) pointing resolution 
(Left) and accuracy (Right). In each case there is little improvement (within 
fluctuations) from beyond 600 hydrophones per cubic kilometre. 

location, the angle out of the plane of the pancake and the incoming neutrino 
energy, is calculated for a range of energies and the energy at which it is a 
minimum should then correspond to the energy of the incoming neutrino. 

CM 

10^ 
10= 



10* 

10 11 12 13 14 15 16 

log^^E (GeV) 

Figure 5.15: Plot of as described by equation l5.15l for one event with suc- 
cessful vertex reconstruction with (solid curve) and without (dashed curve) 
refraction. In each case the energy is reconstructed below the true value. 

As can be seen from the plot in Figure ISTSl the minimisation has yielded a 
result that is below the actual neutrino energy; however, this is easily explained. 
Here the expected energies were calculated merely as function of the distance of 
the resolved vertex without considering the attenuation due to being outside the 
plane of the pancake. So, the expected pulses on those phones outside the pan- 
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cake plane is over-estimated, the true neutrino energy would have to be slightly 
larger in reality to produce these signals. We can see this effect more easily in 
Figure \57[6\ which displays /ogioErecon on the horizontal axis and /ogioEgctuai on 
the vertical axis, for over a thousand events. Superimposed on the plot is the line 
y = X to aid the observation of the systematic shift of the Erecon towards lower 
energy. The mean reconstructed energy is (61 ± 17)% of the actual energy in the 
presence of refraction and (63 ± 12)% without. Naturally one can assume that 
adding the effects of angular attenuation to the calculation will improve the 
reconstruction; however this has not been undertaken in this work. Instead, to 
give an estimate of the potential energy reconstruction performance Erecon has 
been multiplied by a scale factor in order to systematically shift values to higher 
energy, this can be seen plotted in Figure I5.16[ The resulting energy resolution is 
plotted in Figure l5?T7l Up until this point it has been assumed that a constant 25% 
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Figure 5.16: (Left) The distribution of £recon v.s.Eactuai/it is offset from the 
dashed line Erecon = Eactuai due to systematic underestimation of the neu- 
trino energy. (Right) The shifted distribution of Erecon v.s.Ea^tuai compensat- 
ing for the underestimation of the neutrino energy. Both plots include the 
effects of refraction. 



of the energy of the incident neutrino is imparted upon the resulting hadronic cas- 
cade which in turn generates the thermoacoustic pressure wave detected a given 
hydrophone array. As discussed in Section |3^ the proportion of the neutrino en- 
ergy liberated by the cascade can in fact fluctuate from — 100%. Including the 
effects of a non-constant Bjorken-y by sampling from the distribution plotted in 
Figure |3.3(b)| value reduces the effectiveness of the energy reconstruction (by a 
factor orders of magnitude greater than the underestimation error); the resulting 
energy resolution is plotted in Figure \5.17[ 
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Figure 5.17: (Left) Potential energy resolution with c^iogj^E ~ 0.1, assuming 
a constant Bjorken-y of 0.25. (Right) Energy resolution with c^iogjo ^ 
smeared by a non-constant Bjorken-y (dashed) on top of the unsmeared en- 
ergy resolution (solid). Both plots include compensation for systematic un- 
derestimation of the neutrino energy and include the effects of refraction. 



5.9 The Sensitivity Calculation and Predicted Results 

The simulations described thus far determine the ability of a given hypothetical 
array of hydrophones to detect UHE neutrinos across a range of energies. For 
a given energy, it is therefore possible to place an upper limit on the neutrino 
flux, assuming that no events are detected within a certain time window. The 
following calculation is based on the formalism laid out in reference II50II . The 
principal equation is: 

• 0(E) is the limit on the differential flux [dNi, / dEdAdCldt] in units of 
GeV~^cm~'^ sr~^ s~^. 

• Sup is the upper limit on the number of detected events that one places. 
Assuming that no events are seen, then 95% confidence level upper limits 
correspond to setting Sup = 3.0 (or 2.3 for 90% confidence level) in the 
absence of background 1185 II . 

• A(E) is the number of events that we would detect if a unit differential flux 
were incident on our detector at a unique energy E, for time T. It can there- 
fore be written as: 

A(£) = (TyNiE) xpxNAxVxlnx fud^) x T (5.17) 

where: 
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- (7i/]v(E) is the neutrino-nucleon cross section at energy E. This is an 
extrapolated result (see Appendix IA.1|) . 

- |0 is the density of water. 

- Na is Avogadro's number. 

- V is the fiducial volume (i.e. the can) where V = n x (Rcan)^ x Hcan- 

- T is the exposure time in seconds 

- /mc (E) is the fraction of Monte Carlo events generated at energy E that 
are detected (and reconstructed). 

The sensitivity of an array of 1000 hydrophones, distributed at random in a vol- 
ume of one cubic kilometre, is plotted in Figure l5J8l (labelled lyr-35mPa-1000), the 
pressure threshold is 35 mPa and the period of operation one year. Also shown is 
the sensitivity of 1500 km^ array, with 100 hydrophones per cubic kilometre op- 
erating at a pressure threshold of 5 mPa for 5 years (labelled 5yr-5mPa-1500). This 
is an indicator of the size and scale of an array required if one is to be sensitive 
to the cosmogenic neutrino flux (labelled GZK) resulting from the attenuation of 
UHECRs on the 2.7 K CMB photons. The sensitivity curves for this study and 
the other experiments included on the plot represent a limit on the neutrino flux 
based on no neutrinos being detected within the period of operation. Also shown 
are the neutrino fluxes predicted by some theoretical models. 



5.10 Summary 

The details of a simulation dedicated to calculating the ability of hypothetical 
large scale arrays of hydrophones to detect a flux of UHE neutrinos, has been re- 
ported. Through the use of ray tracing algorithms, the effects of refraction have 
been incorporated on an event by event basis on cubic kilometre scales, build- 
ing on previous work ( II80ILII69II ) which neglected or approximated refractive ef- 
fects. The treatment of vertex reconstruction, which in a homogeneous medium 
can be solved analytically with almost perfect accuracy (see Figure I5.3|) , presents 
a non-linear problem, without an analytically determinable solution. Numeri- 
cal computation of refracted vertices has subsequently been discussed, extending 
the work beyond that which incorporates pre-computed look-up tables ( II55ILII67II ) 
which prove to be computationally expensive when looking to simulate arrays of 



90 



5.10. Summary 




5 6 7 8 9 10 11 12 13 14 15 16 

Log[E(GeV)] 



Figure 5.18: Predicted limits on the neutrino flux from this study (A and 
A), RICE (•) and FORTE (■). Also shown are those fluxes predicted by the 
Waxman-Bahcall (dotted). Topological Defect (solid) and GZK (dash-dotted) 
models from Figure IZTSl 

thousands of cubic kilometres. Additionally the dependencies of effective vol- 
ume, vertex resolution and pointing accuracy on hydrophone density have been 
investigated. A full study of array optimisation is required to maximise perfor- 
mance for all three effects in unison. Array performances have been translated 
into limits on the neutrino flux as a means of comparing the sensitivity of dif- 
ferent instrumented volumes. In order to observe and successfully reconstruct a 
so-called GZK neutrino, it has been shown that hydrophone arrays encompassing 
volumes of several thousand cubic kilometres and operating at pressure thresh- 
olds of around 5 mPa are required. 
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Chapter 6 

First Data from Rona Underwater 
Acoustic Range 



6.1 Introduction 



The ACoRNE collaboration has access to the QinetiQ Underwater Acoustic Range 
at Rona in Northwest Scotland. This facility is operated and maintained by Qine- 
tiQ on behalf of the Ministry of Defence (MoD) and is used for the measurement 
of radiated acoustic signatures from Royal Navy ships and submarines in broad 
equivalence to the United States AUTEC array in the Bahamas as used by the 
SAUND collaboration. The array comprises eight wide band hydrophones dis- 
tributed along Inner Sound - a 12 km long strait separating the Inner Hebridean 
islands of Rassay and Rona from the Scottish mainland - with an average depth 
of ^ 230 m. Six of the hydrophones are located at mid-water on cables that are 
anchored at the lower end to the seabed and pulled taut by a buoy tethered at the 
upper end, one is anchored to the floor and the eighth is 30 m from the bottom. 
The array extends over about one and a half kilometres and thus has a sound 
transit time of approximately one second. Since December 2005 data acquisition 
has been underway. The first run constituted 15 days of continuous readout at a 
sampling frequency of 140 kHz and amounted to 2.8 TB of raw information. 
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6.2 Data Acquisition 

6.2.1 Recording the hydrophone signals 

The eight Rona hydrophones are commercially available wide band sensors with 
a flat frequency response up to resonance at 50 kHz. The ACoRNE data acqui- 
sition (DAQ) philosophy is to send all data to shore unfiltered and archive it for 
processing and analysis offline. Every ten minutes each of the Rona hydrophones 
produces a waveform audio format (WAvJ^ file; the eight files are placed in a 
time-stamped directory and written immediately to an 8 TB RAID disk array. In 
September 2006 the DAQ was upgraded to incorporate an LT03 tape drive with 
mechanical loader and barcode reader. The WAV files are compressed using the 
Free Lossless Audio Compression (FLAcj^ codec and written onto 400 GB LT03 
tapes nightly and then removed from the RAID array. Compression via FLAC 
reduces the file size to approximately 45% of its orginal size, with no loss of infor- 
mation. Because of the sensitive nature of the Rona facility remote access to these 
data is impossible; as such tapes are physically mailed to ACoRNE members by 
the operating staff. No filtering whatsoever has been introduced thus far, so prior 
to analysis the raw data remains intact. 

The motivation for such a byte-intensive system is to enable the use of linear- 
phase digital filters. Unless the phase response of a recording system is linear, 
different Fourier components of the signal can lag over different timescales and 
the signal becomes distorted. As such a zero distortion, linear phase DAQ is de- 
sired. Analogue filters that meet these criteria are readily available; however, dig- 
ital signal processing offers a greater degree of versatility and is thus employed 
instead. Digital linear phase filters require a time-reversed signal input in order 
to maintain zero distortion. This cannot be achieved in real time which is why 
unfiltered data are retained. Furthermore there is no limit on the number of times 
the unfiltered data can be reprocessed following any further developments in the 
DAQ system. If a single set of analogue filters were implemented before the data 
were written it would be a one shot mechanism and the raw data would be irre- 
versibly changed. Moreover, it has already been pointed out that several matched 
filters will be required for optimised detection of neutrinos at all angles and dis- 
tances (see Section I4.4|) . Offline analysis of unfiltered data thus allows several 

^http://www.microsoft.com/ Accessed 03/03/2007, http://www.ibm.com/ Accessed 
03/03/2007 

2http://flac.sourceforge.net/ Accessed 03/03/2007 
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passes through multiple matched filters, giving the greatest degree of flexibility 
for matched filter development. 



6.2.2 Reducing the data 

The current DAQ chain is composed of four triggers that are designed to max- 
imise retention of neutrino-like signals. The emphasis of this study is to cat- 
egorise as many types of potential background as possible such that they can 
be readily identified and subtracted from future data, in the hope that what re- 
mains is the signal observed from UHE neutrino induced particle cascades in the 
sea. In order of increasing priority, candidate events were triggered in the initial 
offline data reduction by: pressure (T1=P); the derivative of pressure with re- 
spect to time (Y2=dP / dt); the second derivative of pressure with respect to time 
{T3=d^P/ df-); and a Matched Filter (T4=A1). The raw data were analysed in 10 s 
intervals. For each interval the time stamp for the five most energetic triggers in 
P, dP / dt, d-^P / df- and triggers are taken. 1000 samples (7.14 ms) either side of 
the event are retained for each time stamp and an overlap is said to occur if two 
triggers are within 500 samples. If two triggers of the same type occur within the 
500 sample coincidence window, then the most energetic is taken and the other 
discarded. If triggers of different types occur within the 500 sample window then 
the trigger with the highest priority as described above is retained. 



6.2.3 Rona calibration tones 

Roughly once a day a test of the data transmission between each of the Rona hy- 
drophones and the shore station is performed. A 1 kHz sinusoid is injected elec- 
trically behind the cable-hydrophone connection and the signal is recorded at the 
shore station, to verify the integrity of data transmission along the cable. The cal- 
ibration tones are injected into the system sequentially, that is one hydrophone at 
a time, and are thus easily rejected by a demand for coincidence between events. 
No 1 kHz signals have been observed in this analysis as would be expected since 
a demand for coincidence rejects the sequential signals. 
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6.3 A Novel Analysis of the Reduced Data 

6.3.1 Event discrimination 

A total of 12, 453, 500 signals were retained from the first 15 days of Rona data 
by the data reduction code. An attempt has been made to identify and classify 
different types of event within this data set. In concordance with the large scale 
detector simulation discussed in previous chapters, only those events above a 
threshold pressure of 35 mPa and with a four-fold or greater hit multiplicity were 
considered. In order to discriminate against random noise, Fourier analysis of 
the events was applied such that only those events that showed the same peak 
harmonic in each coincident hydrophone were retained. 

Fourier analysis facilitates identification of the frequency composition of a sig- 
nal. The FFTW algorithm^ was employed to calculate a discrete Fourier transform 
(DFT) of the sampled data. An effective DFT requires a window function iv{n) to 
select a finite number of samples of the raw data x{n). The sampled data can be 
thus expressed by the following equation: 

v{n) = w{n)x(n) (6.1) 

where w(n) is the aforementioned window function. There are many forms of 
window function available, a couple of popular examples being the Harming and 
Blackmann Windows as defined by Equations 16.21 and 16.31 respectively. 

win) = ^[11-cos^] (6.2) 

T-Tzn Atctl 
win) = 0.42 -0.5 cos— — - + 0.08 cos— — - (6.3) 
^ ^ N- 1 N -1 

However, such window functions are best suited to analysis of continuous pe- 
riodic signals and were not implemented. The type of signal one expects from 
a neutrino is much more impulsive. A so-called Rectangular Window is better 
suited to transient signals, this is expressed in Equation |6.4| below: 

w{n) = 1 for 0<n<N-l 

w{n) = elsewhere. (6.4) 

In general a DFT requires that the number of samples is an integer power of 
two. In this analysis, an event was viewed through a window of N = 2"^ = 128 



3http://www.fftw.org/ Accessed 13/03/2007 
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samples, corresponding to width of 128/140000 = 914.28 las; this is the sample 
period. Remember, a typical bipolar signal from a neutrino with a dominant 
frequency of 10 kHz produces a pulse ~ 100 las wide and thus fits comfortably 
within this window. A smaller window could have been chosen, to better iso- 
late events, but at a cost of reducing the number of harmonics identifiable in 
the DFT. There are N/2 = 64 fundamental harmonics of the sampling period 
in a DFT of N = 128 samples of data; the nth harmonic is at n x 140000/128 = 
n X 1093.75 Hz. A window of N = 128 samples was chosen because it gave a 
good compromise between signal isolation and harmonic identification. 

Identification of events showing the same dominant frequency was achieved 
by requiring that all coincident events showed a maximum in the same bin (i.e 
harmonic) of the DFT histogram. This is illustrated clearly by the identifica- 
tion of sinusoidally oscillating tones present in the data, as plotted in Figure 16.11 
Discrimination of sinusoidal signals is aided by the fact that they show distinct. 




Figure 6.1: Example of a 10 kHz oscillating signal. The left hand display 
shows the pressure signal in mPa and the right hand plot shows the mag- 
nitude of the DFT. In the pressure signal the x-axis corresponds to sample 
number, there are 128 samples in the event window producing 64 harmonics 
in the DFT. A clear peak in the tenth harmonic of the DFT is indicative of a 
single frequency sinusoidal composition at ~ 10 kHz. 

isolated peaks in the DFT. If a sinusoidally oscillating signal has a period that is 
an integer number of sampling periods, a peak occupying a single bin in the DFT 
would be observed. However, in the case plotted in Figure [Ol the signal period 
is 100 las and the sampling period is 128/140000 = 914.29 las. This mismatch is 
the source of the tails either side of the peak in the tenth harmonic in the DFT; the 
effect is known as leakage. 

6.3.2 Categorisation of events 

Of the 12, 453, 500 reduced events, 3526 corresponded to multiplicities of four or 
greater above the 35 mPa threshold with an equivalent frequency composition as 
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identified in the DFT. The number of triggers per hydrophone and the trigger 
frequencies are plotted in Figure I6.2[ 




Figure 6.2: The number of triggers before (unfilled) and after (filled) cuts on 
each hydrophone at Rona, as retained by the data reduction code (Left). The 
number of instances if each trigger (Right). 

Oscillating signals 

One hundred and twenty-nine of the 3526 events that passed the selection criteria 
were signals of a sinusoidally oscillating nature. These events, whilst originally 
being classified as calibration tones, did not occur at occasions when a Rona cal- 
ibration tone was scheduled to be injected into the system, as recorded in the 
recorder logbook entries. Nor were the frequencies consistent with those injected 
by the Rona operators. An example of such a signal is plotted in Figure 16. 3[ The 
source of these signals is currently unknown but it is suggested that they repre- 
sent an anthropogenic background as opposed to bio-acoustical or environmental 
noise. All of the 10 kHz oscillating signals identified in this analysis occurred on 
14/12/2005 between 11:24 and 14:14. Two higher frequency oscillating signals 
were recorded with a five-fold coincidence at 10:24 on 14/12/2005 at 10:24, one 
high frequency signal on 16/12 at 04:44 and one further high frequency oscillating 
signal on 24/12 at 07:33. An example of one of the high frequency hydrophone 
traces is plotted in Figure 16.41 the remainder of the oscillating signals being at 
10 kHz. 

Transient events 

The remaining 3397 events were further subcategorised in terms of their fre- 
quency composition. Firstly the DFT of an event was analysed and classified 



98 



6.3. A Novel Analysis of the Reduced Data 



as one of three types: smooth, mixed, or peaky as according to the criteria in Table 
16.11 The threshold for peakiness of 2.5 is somewhat arbitrary but was chosen be- 



Peakiness 


Mean Number of Peaks in DFT 


smooth 


Speaks < 2.5 


mixed 


Speaks 2.5 


peaky 


Npeaks > 2.5 



Table 6.1: Events were categorised in terms of their Peakiness, a variable 
determined by the number of peaks in the DFT. 

cause it gave a good separation between events that showed a continuous DFT 
spectrum and those that showed a relative separation between peaks in the DFT. 
An example of two such events is plotted in Figures [631 and [63! 
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Figure 6.3: A 10 kHz event coincident on six hydrophones showing a clear 
sine wave component. Traces are Pressure (Pa) (Left) and DFT magnitude 
(Right). Sample number runs along the x-axis. 




Figure 6.4: Example of a high frequency signal at Rona. The peak of the DFT 
is in the 46th harmonic, corresponding to a frequency of ~ 50 kHz. 
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Observing the DFTs of the reduced Rona events, it is clear to see that there are 
two extremes in the acoustic data. The first is a periodic, oscillating signal such 
as the 10 kHz sinusoid in Figure [63l This type of signal produces a distinct single 
peak in the DFT The second is a type of signal evident in Figure [631 This appears 
to be an impulsive event, as indicated by the continuous DFlQ. The pulse shape 
of impulsive events is dominated by the impulse response of the hydrophone 
system whereas continuous oscillating signals show little contribution from the 
detector elements. Studying the shape of the DFT, i.e. their peakiness, thus gives 
an insight on the nature and origin of these events. 

Analysis of the remaining 3397 events that did not categorise as sinusoidally 
oscillating signals yielded 2329 smooth, 44 mixed and 969 ipeaky signals. These events 
were further sub-categorised in terms of the magnitude of their frequency com- 
ponents according to the criteria in Table 16.21 No events satisfied the criteria for 



Frequency Range (kHz) 



low ( < f < 22) 

medium (22 < f < 44) 
high (44 < f < 70) 



Mean DFT Harmonic 



< DFT < 20 



20 < DFT < 40 



40 < DFT < 64 



Table 6.2: Rona events were sub-categorised according to the mean value of 
the DFT histogram. 

classification as high frequency. Only 55 events classified as medium frequency 
events thus indicating a clear dominance of events below ^ 22 kHz. 



6.3.3 Time dependence of events 

Investigation into the temporal distribution of Rona events passing pressure and 
coincidence cuts suggests that activity is intermittent. There can be periods of 
sustained activity, lasting for several hours, showing events in every consecutive 
ten minute WAV file; and, there can be periods where a high number of events 
occur within a very short space of time, occupying only a single ten minute WAV 
file. The temporal variation in activity is best illustrated in a plot of event fre- 
quency as a function of time. Figure [bTI shows this. 

*the DFT of an impulse is a constant: DFT{5{n)} = 1 
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(a) 2D temporal distribution of Rona events 
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(b) 3D temporal distribution of Rona events 

Figure 6.7: 2D and 3D plots of the temporal distribution of Rona events. Ac- 
tivity appears to be constrained to bursts of varying length and intensity. A 
large number of events can occur in a short space of time, producing isolated 
peaks such as the one at 21:00 on the 10th; or, events can be produced during 
a burst sustained over several hours such as the period betv\^een 08.00 and 
16.00 on the 14th. 
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6.3.4 Recurring event types 

Certain types of event, with distinct characteristics, are seen in the data on a recur- 
ring basis. An attempt has been made to categorise such events. Six distinct sig- 
nal shapes have been identified. What follows is a brief description of each event 
type, ending with a table summarising the dates on which each of the events 
occur in the Rona data. For each type of event a single hydrophone signal inclu- 
sive of DFT is plotted, pressure in Pa runs along the y-axes and sample number 
along the x-axes. The full event displays showing coincidences can be found in 
Appendix El 



Short impulsive 
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Figure 6.8: A short-lived impulsive event. 



As already introduced in Section [6.3.2l there is a category of event that appears 
to display an impulsive component. It is categorised by a short duration multi- 
polar waveform and a continuous DFT. It is however indeterminable from the 
DFT alone whether the signal is actually the response of a hydrophone to an 
impulse or if it really is a bipolar signal with a ~ 20 kHz peak frequency as one 
might expect from an UHE hadronic cascade. 



Low frequency oscillator 





Figure 6.9: A low frequency oscillating event. 



An event categorised by a low frequency (< 4 kHz) oscillating waveform ap- 
pears consistently throughout the first fifteen days of Rona data. The pressure 
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wave of these signals can often appear as in Figure IE.2I or it can be somewhat 
distorted. Nevertheless this kind of signal is present in almost every day's worth 
of data (see Table [63)) and thus constitutes a significant background. 
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Figure 6.10: A 10 kHz sinusoidally oscillating event. 



Again, this signal has previously been introduced in Section 16.3.11 At first it 
was thought that these 10 kHz oscillating signals were in fact calibration charges 
injected into the Rona DAQ to verify the integrity of transmission along the eight 
hydrophone to shore cables. However, this was later ruled out because they do 
not occur at times when the calibration tones were injected nor are they at the 
right frequency. Nonetheless, it is suggested that these signals are anthropogenic 
in origin. 

High frequency oscillator 




Figure 6.11: A high frequency event. 



Higher frequency oscillating signals are also present in the Rona data but with 
a much lower abundance than the low frequency oscillators and the 10 kHz sinu- 
soids. Typically they occur at frequencies around 50 kHz. 

Ringer 

Some events in the Rona data occur periodically that demonstrate a decaying and 
oscillating pressure signal. There are two possible reasons for the appearance of 
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Figure 6.12: A 'ringing' event. 



such a signal: firstly, there is a decaying oscillating pressure wave in the water; 
secondly, due to a non-linearity in the phase response of the hydrophone system 
an asymmetric time delay of different Fourier components has caused the signal 
to ring. It is impossible to know for certain which of the two explanations is cor- 
rect for any given case so these events have been arbitrarily labelled as "ringers". 

Bipolar 




Figure 6.13: A 'bipolar ' event. 



Finally, this analysis has identified a class of event that is approximately bipo- 
lar in shape with a frequency close to that which we would expect from a neutrino 
induced particle cascade at some far-field distance from the Rona array. The sig- 
nal is certainly not flat either side of the bipolar part of the pulse, and may indeed 
contain an oscillating component but one can imagine that if a smaller window 
had been chosen a clearer bipolar signal would be observed. 

Clearly there are noise sources around the Rona array that can portray them- 
selves as originating from UHE hadronic cascades. Further, more detailed work, 
beyond the scope of this somewhat qualitative analysis, is required to fully un- 
derstand such backgrounds. 

Occurrences of identified signal types 

The identified signal types as described above recurred on several days through- 
out the data analysis. Table [63] summarises this information. 
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Event Type 


Recurrence 


Short impulsive 


11th, 12th, 14th, 20th 


Low frequency oscillator 


9th-12th, 14th-17th, 19th-21st, 23rd 


10 kHz sinusoidal 


14th 


High frequency 


14th, 16th, 24th 


Ringer 


12th, 21st, 24th 


Bipolar 


17th, 19th-22nd, 24th 



Table 6.3: Recurrence of events at Rona. 



6.4 Summary 

The DAQ as was in operation between 09/12 and 24/12/2005 at the Rona hy- 
drophone array in Scotland has been described. Additionally a novel analysis of 
the raw 12, 453, 500 events triggered by the ACoRNE data reduction has been re- 
ported. 3526 passed offline cuts demanding that: there is a coincidence between 
four or more hydrophones; the peak pressure is above 35mPa; and, the maxi- 
mum peak in the DFT of each coincident signal is at the same frequency. The re- 
quirement for four-fold or greater coincidence above a threshold of 35 mPa puts 
the analysis in concordance with the sensitivity simulation described in previous 
chapters. Table 16.41 summarises the scale of coincidence amongst the analysed 
data. A total of 1059 events showed a five-fold coincidence or greater. In Section 
14.41 it was stated that for a mean noise level of 35.9 dB only one false trigger in 
ten years, due to noise, with a five-fold coincidence should be observed. By in- 
tegrating the energy of the pulses triggering the ACoRNE DAQ in this 15 days 
worth of data, a noise level of between 90 100 dB is observed, which demon- 
strates a higher level of noise activity than the location considered in the ACoRNE 
proposal. Furthermore it is suggested that there is an anthropogenic component 
of the Rona events that raises the noise level above that predicted by sea state 
alone. The noise at Rona provides an abundance of "signals", the observation of 
which enables development of the data acquisition system. The more developed 
DAQ can then, in theory, be implemented elsewhere at a warmer, deeper and 
ultimately quieter site more suited to the acoustic detection of UHE neutrinos. 



Level 


4-fold 


5-fold 


6-fold 


7-fold 


8-fold 


Frequency 


2468 


845 


182 


32 






Table 6.4: Coincidences at Rona. 
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Chapter 7 

Discussion and Outlook 



7.1 The Motivation for Neutrino Astronomy 

The motivation for Neutrino Astronomy is broad and compelling. As was dis- 
cussed in the introduction the eventual first detection of neutrinos lagged their 
initial postulation by some twenty-three years. They are an essential component 
of the Standard Model of physics and are the second-most abundant known par- 
ticle in the Universe yet their properties are still shrouded in mystery. One thing 
is for sure: the UHE pp and pj interactions that follow the production of the 
highest energy cosmic rays will carry with them an associated flux of UHE neu- 
trinos (Section II. 3. Due to their infinitesimal mass and lack of electric charge 



they can offer a unique perspective of the Cosmos. What is more, the scope for 
so-called "new physics" offered by Neutrino Astronomy is abundant, whether it 
is from the decay of exotics such as topological defects or relic big bang particles, 
or the annihilation of dark matter, neutrinos can offer a window into the frontier 
territories of quantum scale gravity, extra dimensions, super-symmetry and dark 
matter (Section [L3]2]). 



7.2 The Current Scale of Development 

The Neutrino Astronomy effort is a global endeavour extending across the sur- 
face of the Earth from pole-to-pole and above-and-below from high in the atmo- 
sphere to the depths of the sea (Sections 12.31 12.41 and |2.5|| . It is because of the 
pioneering efforts of those such as the ill-fated DUMAND project and the long 
established lake Baikal experiment that today's neutrino telescopes can look to 
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the future and envisage a technology on the cubic kilometre scale and beyond. 
Despite such grand aspirations the field as a whole remains in its infancy when 
compared, say, to the status of radio astronomy and to a lesser extent UHECR and 
TeV gamma astronomy One hopes that in the years to come the fruits borne by 
this labour can match that which is emerging from experiments such as AUGER 
and HESS in the present day. 

The acoustic detection of UHE neutrinos is undoubtedly the least well devel- 
oped technique. However it is arguably the most elegant in its simplicity. The 
direct detection of the thermoacoustic emission resulting from the interaction of 
UHE neutrinos at the Earth, with an off the shelf technology such as commercially 
available wide band hydrophones, is unrivalled in its potential to offer huge ef- 
fective volumes without the necessity of complicated detector design. 



7.3 Simulation of Acoustic Arrays 
7.3.1 Neutrino interactions 

The ACoRNE simulation work has continued to be developed following the for- 
mation of the collaboration. From calculations of the initial interaction of a neu- 
trino at the Earth to the reconstruction of the thermoacoustic emission by cubic 
kilometre arrays of hydrophones a vast number of processes require considera- 
tion. Understanding the neutrino-nucleon cross-section at EeV energies depends 
on some combination of parameterised empirical data, extrapolated to higher 
energies, and theoretically motivated approximations. Man-made particle accel- 
erators will likely never rival those of nature in their ability to accelerate parti- 
cles to the highest energies. Presently it is from CR air shower experiments that 
we gain the most insight into particle interaction cross-sections at the TeV scale 
and beyond. Eventually neutrino telescopes, it is hoped, will provide the same 
quality of information as experiments such as HI-RES, KASCADE, AUGER and 
those alike, allowing the existing models that are relied upon in programs such 
as ANIS, PYTHIA and CORSIKA to be tuned to a finer degree of accuracy. For 
now we accept that the level of uncertainty in our particle interaction simulations 
is at worst within the level of individual event fluctuations (Sections 13.3.21 and 
Appendix |B]). 
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7.3.2 Formation and propagation of the signal 

The formation of the bipolar acoustic signal characteristic of an UHE neutrino 
induced particle cascade relies exclusively on knowledge of the thermoacoustic 
energy deposited by the hadronic cascade. Efforts both experimental and theoret- 
ical have been made to best understand this mechanism. In the simulation work 
contained in this thesis, parameterisations of the integrated thermal energy re- 
sulting from the Geant4 and CORSIKA programs at various energies have been 
utilised to these ends (Section |3.4[) . Furthermore, the far-field approximation has 
been invoked to parameterise the peak pulse pressure, which has been shown to 
vary in direct proportion to the energy of the cascade. 

The propagation of acoustic signals in the sea is as rich and diverse a field as 
Neutrino Astronomy. To attempt to incorporate something with such a magni- 
tude of complexity into simple parameterisations of absorption, attenuation and 
propagation (Sections 14.21 and |43l) is perhaps something of an injustice. Neverthe- 
less progress will continue to be made upon the first generation codes that incor- 
porated a constant velocity of sound towards more complicated ray-traced meth- 
ods as described in this work, complex absorption models and beyond. Knowl- 
edge of the sound velocity profile surrounding a detector is critical to our under- 
standing of its performance. Indeed we are fortunate that the deep sea presents 
something of an idealised situation in that the SVP is essentially linear, yet it re- 
mains to be seen exactly how accurate our understanding of it need be. What is 
more the boundary conditions at the sea surface and sea bed that have been ne- 
glected in this work present in themselves a broad subject for study. If the bipolar 
acoustic signal can undergo surface or sea bed reflections without experiencing 
too much distortion then the effective volume for a given array may extend well 
beyond the cutoff of 48 km observed in this work. 

7.3.3 Noise and filtering 

The effects of noise and filtering again encompass a vast amount of research that 
is parameterised into the simple models utilised in the large scale detector simula- 
tion described in this thesis (Section |4.4[) . What is being learnt from experimental 
projects such as SAUND in the Bahamas, the acoustic branches of the ANTARES 
and NEMO collaborations at their respective sites, and from the ACoRNE project 
at Rona (Chapter ^ is that the backgrounds presented to the acoustic neutrino 
astronomer are numerous and multifaceted, presenting in themselves events that 
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emulate the type of signal one may expect from an UHE neutrino induced particle 
cascade. Understanding of thresholds continues to develop as data are gathered 
from arrays such as those operated by the aforementioned groups and so does 
the continuing development of matched filters. 

7.3.4 The large scale detector simulation 

The reconstruction algorithms developed as part of this work represent the first 
iteration in what will undoubtedly mature into a vastly more complex field of 
study. Sensitivity predictions for acoustic arrays at present represent the work of 
a handful of people working broadly in isolation. In addition to the development 
of the computerised models, described above, that effectively sit as a front-end to 
the large scale simulation there is scope for further refinement of the techniques, 
already in use in this study, with respect to vertexing, pointing and calorimetry 
(Sections I5.4[ 15.61 and I5.8 [) . Better development of the minimisation procedure has 
already been alluded to (Section |5A2l), additionally optimising the reconstruction 
of the acoustic pancake remains an interesting and open matter for study. One can 
envisage much improved performance by fully exploiting the unique pancake 
geometry. 

The nature of this work, as it was conceived in the beginning, was to consider 
typically one thousand hydrophones in a volume of one cubic kilometre. We 
have seen that there are different factors for which array density can be optimised 
such as the effective volume and the reconstruction performance (Sections 15.51 and 
15. 6[) . Naturally if one were to begin a proposal to design and build a dedicated 
acoustic array for the purpose of UHE neutrino astronomy extra considerations 
such as cost and deployability come into play; the randomly distributed arrays 
investigated here may practically be very difficult to manufacture. Clearly there 
is a great deal of research potential for proper optimisation of the geometry of an 
acoustic array. 

With respect to the sensitivity predictions made here (Section 15. 9|) , there is 
broad agreement with the work published in reference [6911, (which recursively 
agrees with predictions from the IceCube coUaborationjH In loose terms, as cur- 
rent predictions stand, one requires a pressure threshold around 5 — 10 mPa op- 
erating across an array of over a thousand cubic kilometres with at least one hun- 
dred hydrophones per cubic kilometre if the proposed flux of GZK neutrinos is 

^see the presentations and proceedings of the ARENA 2005/06 workshops 1861, 1871 
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to be observed. Given further study and development of array optimisation and 
reconstruction algorithms, this may change. 

7.4 Predictions for Future Work 

Within the field of acoustic neutrino astronomy there is now some idea of the 
scales (thousands of cubic kilometres) to be considered if one wants to guarantee 
observation of those neutrinos resulting from the attenuation of UHECRs by the 
cosmic microwave background radiation. Undoubtedly that level of commitment 
is years beyond today's level of activity. Currently, acoustic projects exist by util- 
ising existing infrastructures, be it through military hardware such as the AUTEC 
and Rona arrays used by SAUND and ACoRNE respectively, or by piggybacking 
the deployment of existing sea or ice based detectors such as ANTARES, NEMO 
and IceCube. This latter strategy paves the way for hybrid detection. Because 
water doesn't facilitate radio detection and there is no overlap in energy for si- 
multaneous optical and acoustic detection of a given neutrino (see Figure I2.15|) , 
ice is surely the best medium for hybrid detection. This is further enhanced by 
its higher density and longer light absorption lengths than for water, coupled 
with the potential for shear wave acoustics thanks to the lattice-like structure of 
ice. Within the IceCube collaboration work is underway to extend to hybrid op- 
tical/radio/acoustic detectiorM Any detector offering simultaneous registration 
of the same event via two independent means has intrinsically less systematic un- 
certainty. One as yet unanswered question that remains for the acoustic neutrino 
astronomer is that of positively identifying a signal as being that of a neutrino 
and completely ruling out any other source. Naturally a coincident radio based 
detection would remove a great deal of uncertainty. 

As the ANTARES project nears full deployment the boreal cubic kilometre 
scale neutrino telescopes begin their design phases. Acoustic instrumentation is 
already a part of the ANTARES design scheme and is almost certainly going to 
feature in the next generation telescopes in the Mediterranean. If this existing 
infrastructure can be utilised in such a way as to provide readout for acoustic 
detection then the first cubic kilometre, sea-based acoustic neutrino telescopes 
may not be far away thus completing complementary neutrino telescopy in both 
the northern and southern hemispheres, a requisite for full sky coverage. The 

^again, see for example the proceedings of the ARENA 2005/06 workshops Il86l , ll87ll 
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neutrino-light in which the Earth bathes will soon illuminate the world's largest 
active particle detectors and enlighten physicists to the mysteries foretold by na- 
ture's most enigmatic cosmic messenger. 
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Appendix A 
Empirical Data 



A.l Neutrino-Nucleon Cross Section 




Figure A.l: The total neutrino-nucleon cross section and its charged current 
and neutral current contributions. Taken from II88II . The cross section for 
neutrino interactions at UHE was extrapolated from this plot and translated 
into the following code: 



// logl0(sigma/cin-2) = -32.60 + (0.362*(loglO(E/GeV)-7)) 
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// INPUT : energy in eV 

// OUTPUT : cross section in cni"2 

double egev = energy/lE09; 

double loglOcs = -32.60 + (0.362*(logl0(egev)-7.0)) ; 
return pow(10 , loglOcs) ; 



A.2 Knudsen Noise Curves 
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Figure A.2: Typical Knudsen noise curves (SSO, SSI, SS3) at the RONA site 
off the west Scottish coast used by the ACoRNE collaboration. These curves 
parameterise the ambient noise levels as a function of Sea States as defined 
by the World Meteorological Organisation. Also included is the measured 
noise level at Rona (Refer to legend). 
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Comparison of Hadronic Showers 
from GEANT4 and CORSIKA 

The longitudinal and radial profiles of thermal energy deposition from hadronic 
cascades produced by two programs: GEANT4 and CORSIKA are compared. The 
CORSIKA program, designed for the production of UHECR induced air show- 
ers, in the upper atmosphere of Earth, has been modified to produce hadronic 
cascades resulting from UHE neutrinos in a seawater medium |[72||. 

In each case the average of 100 showers is plotted to smooth out fluctuations 
on an event by event basis. The agreement between the shower shapes from 
GEANT4 and CORSIKA is within the 20% 30% fluctuations seen between 
showers and the ~ 30% uncertaintity between choice of theoretical models in 
each program. 

The modified version of CORSIKA has been developed so thermal energy 
densities for hadronic cascades initiated by primaries beyond the 100 TeV limit 
inherent to GEANT4 can be produced. Hence the extrapolated pressure to energy 
parameterisation need not be relied upon. A study of shower shapes and peak 
pulse pressure resulting from CORSIKA showers has however validated the as- 
sumption that pressure scales linearly with shower energy. 
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Figure B.2: Geant4 and CORSIKA longitudinal energy deposition at 10^ eV. 
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Figure B.3: Geant4 and CORSIKA radial energy deposition at 10'' eV for 
various depths along the shower axis. 
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Figure B.4: Geant4 and CORSIKA radial energy deposition at 10^ eV for 
various depths along the shower axis. 
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Appendix C 

Formalism of Ray Tracing Algorithm 



We start with the differential equations for ray paths in inhomogeneous (meaning 
the velocity of sound is not constant) media as given in 1182 II and proceed through 
the derivation therein: 

ds\ dsJ dx 
ds ^ ds ) dy 

— (n—) = ^ (C3) 
dsKdsJdz 

where index of refraction n = Cn/c and c„ is the velocity when the angle 9 of the 
ray with respect to the horizontal is 0. If the velocity c is a function of z only then 
n = n{z) and Equations IC.ll through IC.31 are reduced to: 

dx 

n— = constant 

ds 

dy 

n-f^ = constant 
ds 

d fdz\ _ dn 
ds\ dsJ dz 

We can arbitrarily collapse into 2 dimensions by observing that the first two equa- 
tions imply the ray path lies in a plane orthogonal to the xy-plane. From Figure 
IC.ll we see: 

dx 

— = cos 6 



ds 

dz 

ds 



sine (C.5) 



121 



Appendix C. Formalism of Ray Tracing Algorithm 



X 




z 

Figure C.l: Geometry of a single curved ray in a medium with varying 
sound velocity. 



Taking the first equation of Equation IC .41 and substituting the first equation from 
Equation IC.5I 

^ = 1 = K (C.6) 

C Cfi 

which is Snell's Law. K is called the ray parameter which, for the geometry in 
Figure ICT] can be conveniently written as: 

X=^ = ^ (C.7) 

Co Cl 

where Co,Oq and Ci,6i are the sound velocity and angle at (xq/Zq) and (xi,zi) 
respectively. Substituting the second equation from Equation IC .51 into the second 
equation in Equation lC.4l we have: 

dn d / ■ 

Tz = 

d sin 6 dn 

= n — hsm^— 

ds ds 

do dn 
= ncos0— + sm0— 
ds ds 

AO . dn dz 
= ncosfc*— + smfc'— — 
ds dz ds 

^dO . 2^dn 
= ncosO— + smO— 
ds dz 
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which can be solved for dd/ ds, 

do cos 9 dn 
ds n dz 

cos 6 dc 



c dz 

- < <^-^) 

This yields the important result that the curvature of the ray d6/ds is directly 
proportional to the velocity gradient dc/ dz. From the geometry in Figure ICT] we 
can see that 

^ = cote (C.9) 

dz 

such that horizontal range of the ray is given by: 

dx = cot 6dz 

xi-xq = r\otedz (C.IO) 

Jzq 

The travel time along the ray is given by: 

ds 

dt = — 
c 



h - to 



2i ds 



2i dy 

h-to = I (C.ll) 



Izo csm9 

Similarly, from the geometry in Figure IC.ll the ray path length is given by: 

dz 



ds 



sin 9 



/■^i dz 

si-so = / ^ (C.12) 
Jzo sm9 



In the case of a linear sound velocity gradient: 

C = Cni + gZ 

dc 



123 



Appendix C. Formalism of Ray Tracing Algorithm 



where the second term is the derivative with respect to z of the first. We can now 
substitute the differential form of E quation IC . 1 31 into Equation IC .81 

fs - <^-"' 
do / ds is the curvature of the ray, and for a linear sound velocity gradient is con- 
stant. Hence the ray trajectory follows the arc of a circle. The radius of curvature 
R is given by: 

Equations IC.lOl IC.llI and IC.12I for horizontal range, travel time and path length 
respectively can be expressed in terms of 6. Starting with Equation lC.lOl we have: 

Xl — xq = / cot 6dz 

Jzq 

^1 cos 9 
Zo sin0 

^1 cos Odd 



Kc' 

where we have used dz/ sin 6 = —dO/Kc' obtained by differentiating Snell's Law: 

cos 9 = Kc 



dcos9 
dz 

sin 9 d9 



Kc' 



-Kc' 



dz 

dz d9 



sm9 Kc' 

The integral now becomes: 



Xl — Xq 



^1 COS 9 d9 

do Kg 



-^[sm9f' 
Kg^ J^o 

1 

-;— (sin 9o — sin 9-\ ' 

Co ( sin 00 — sin 9\ 
cos 00 
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Similarly, starting with Equation IC. Ill 



h — to 



dz 



zq csind 

de 



00 g<^osd 



1 , I sec 9i + tan 9i 

- — In < — 

g I sec do + tan do 

1, rtan(^ 

- In < -. 

S I tan (^2^ 

1 f cos0o(l + sin^i) 
g I cos 01 (1 + sin^o) 



Finally from Equation IC. 121 



Si - So 



^1 dz 
zq sin6 

1 de 

Oo Kg 



—Oo - 01 
Kg 



° (^0 - 0\] 



gcos 00 



One can formulate an equation for the change in depth of the ray directly from 
Equation IC.13l by using Snell's Law as given in Equation lC.71 From: 



C Cfii 

z = 



8 
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One can write: 

C\ - Cm\ I Cq- Cm 



Z\ - Zq 



Cl - Cq 
g 



1 / COS 01 

- Co — - Co 

g \ COS ^0 



(cos 9i — COS ^o) 



g COS 00 

In summary, when the sound velocity gradient is constant, the horizontal range, 
travel time, ray path length and change in depth are given by: 

Co / sin^o -sin0i\ 

xi-xo = — (C.16) 

g \ COS0O J 

, , 1, /cos0o(l + sin0i)l 

t^-to = — In <^ ) . ' } (C.17) 

g |^cos0i(l + sm^o) J 

si-so = ^V(^o-^i) (C.18) 
g cos 00 

zi-zo = — ^(cos0i-cos0o) (C.19) 
g cos 00 

Equations lC.17[|C.18llC.19l and lC.19l form the basis for the ray tracing algorithm. 
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Neutrino Source Data 



The neutrino spectrum Ey, source distribution {p,6,(p)v and interaction coordi- 
nates (r,z, ^)int as generated in the large scale detector simulation (Chapter [S]) is 
summarised by the following plots (with energy in GeV, angles in radians and 
distances in metres): 
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Recurring Event Types at Rona 

The coincident signals for each event introduced in Section 16.3.41 are as 
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Figure E.2: A low frequency oscillating event. 
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Figure E.3: A 10 kHz sinusoidally oscillating event. 
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Figure E.4: A high frequency event. 
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